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1.  IntroductUMi,  Scope  and  Objectives 


This  report,  titled  "Probability  Based  Ship  Design  Procedures  -  a  Demonstration",  is 
the  second  in  the  series  of  projects  undertaken  by  the  Ship  Structure  Committee  in  the 
thrust  area  of  reliability  based  ship  design.  The  fust  was  the  development  of  a 
comprehensive  primer  to  structural  reliability  theory  as  implied  to  ships  and  marine 
structures.  Ref.  6.  The  work  in  this  project  assumes  that  the  reader  is  familiar  with  the 
various  concepts  and  applications  discussed  in  Ref.  6,  "An  Introduction  to  Structural 
Reliability  Theory",  SSC  Report  351. 

The  immediate  objective  of  this  project  is  to  provide  a  demonstration  of  the  use  of 
probability-based  ship  design  methods  and  to  compare  the  results  with  traditional  design 
methods.  Based  on  the  results  of  the  demonstradon,  the  following  conclusions  and 
infoimadon  are  provided: 

1.  The  benefits  and  drawbacks  of  the  use  of  probability-based  design  methods  compared 
to  the  tradidonal  methods 

2.  The  additional  information  necessary  to  conduct  probabUity-based  ship  designs 

3.  A  summary  of  the  proposed  probability-based  method  showing  how  it  can  be  applied 
to  generate  new  designs  of  uniform  safety  and  how  it  can  be  used  to  assess  the  safety 
of  an  existing  design 

4.  A  discussion  of  the  current  and  future  SSC  projects  in  reliability  and  loads. 


Two  basic  demonstrations  are  provided  in  this  report  (Part  1  and  Part  2)  together  with 
reliability  process  definitions  (Part  3).  These  are  summarized  as  follows: 

1.  Probability-based  design  procedure  -  code  calibration: 

The  objective  of  this  part  is  to  provide  an  illustration  of  how  probability-based 
methods  can  be  used  to  develop  and  calibrate  a  code  (or  design  criteria)  in  order  to 
produce  designs  with  uniform  safety  over  a  wide  range  of  the  basic  parameters  involved 
in  the  design.  For  this  purpose,  ABS  primary  hull  girder  longitudinal  strength  criterion  is 
considered.  A  formulation  for  the  minimum  required  section  modulus  that  satisfies  this 
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requirement  (uniform  safety)  is  developed.  A  demonstration  is  made  of  how  partial 
safety  factors  are  determined,  calibrated,  and  used  in  new  designs  that  have  uniform 
safety. 

2.  Probability- based  ship  safety  analysis: 

The  objective  of  this  part  is  to  provide  an  illustration  of  how  to  apply  state-of-the-art 
reliability  techniques  in  order  to  determine  the  safety  level  of  an  existing  ship  or  an 
existing  de'^.gn,  i.e.,  to  develop  the  ship  safety  indices  taking  into  consideration  the 
uncertainties  associated  with  the  environment,  loads,  matt.nals  and  analytical  nuxlels. 
For  this  purpose  a  tanker  was  selected  in  consultation  with  the  Project  Technical 
Committee  (PTCT)  for  use  in  an  example  to  illustrate  the  safety  assessment  procedure. 
Several  limit  states  were  formulated,  namely  ultimate,  serviceability,  and  fatigue  limit 
states,  and  applied  to  the  tanker.  The  loads  corresponding  to  these  limit  states  were 
developed  and  a  safety  index  was  calculated  for  each  limit  state  using  both  first  and 
second  order  reliability  methods. 


3.  Structural  reliability  process  definitions: 

An  extension  of  the  work  of  this  project  (SR- 1330)  was  approved  by  the  PTC. 

The  additional  work  is  described  in  the  following  tasks: 

(a)  Definition  of  terminology  associated  with  structural  reliability  of  ships  and  offshore 
structures.  This  includes  terminology  related  to  loads,  strength  and  structural 
reliability. 

(b)  Identification  and  description  of  appropriato  ultimate  limit  states  associated  with 
lifetime  extreme  design  loads.  These  include  global  (hull  girder)  initial  yield,  fully 
plastic  and  collapse  limit  states,  and  local  ones  related  to  column,  beam/column  and 
torsional/flexural  buckling  of  longitudinals,  and  grillage  buckling  of  longitudinals 
together  with  transverse  beams. 

(c)  Identification  and  description  of .  rviceability  limit  states  associated  with  plate 
buckling  and  fatigue. 

(d)  A  review  of  probabilistic  extrapolation  techniques  for  lifetime  extreme  loads. . 
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A  NOTE  ON  NOTATION 


A  distinction  needs  to  be  made  between  nmdom  variables  and  their  characteristic  or 
nominal  values,  although  this  may  often  be  evident  from  the  context  In  this  report 
where  necessary,  random  variables  are  denoted  with  a  'tilde'  on  the  top,  e.g.  Cy  is  a 
random  variable,  while  Oy  is  a  nominal  or  characteristic  value. 
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2.  Preliminary  Assessment  of  Reliability  Levels  Implied  in  ABS  Rules 

As  a  demonstration  of  a  probability-based  calibradon  procedure  of  a  code,  the  safety 
level  implied  in  ABS  Rules  for  hull  girder  longitudinal  strength  is  determined  by 
calculating  the  reliability  indices  (P's)  for  300  ships  designed  according  to  the  Rules. 
The  range  of  safety  (Pnmge)  calculated  as  the  difference  between  the  largest 

and  smallest  safety  indices  of  all  the  designs  considered.  An  average  safety  index  (Pay) 
was  also  calculated.  The  objective  of  the  calibration  process  is  to  determine  partial  safety 
factors  to  be  used  in  a  modified  formulation  for  longitudinal  strength  such  that  the 
resulting  safety  level  of  all  designs  is  approximately  constant  with  a  value  equal  to  Pav 
and  such  that  the  resuldng  safety  range  (Piange)  aniong  the  new  designs  is  minimum. 
The  details  of  the  calibration  process  is  illustrated  in  the  following  sections. 

2.1  Limit  State  Formulation 

The  section  modulus  requirements  for  a  ship  according  to  ABS  Rules  is  based  on  a 
permissible  stress  which  is  based  on  the  yield  strength  of  the  material.  For  this  reason, 
only  the  inidal  yield  limit  state  will  be  formulated  which  is  similar  to  ABS  minimum 
secdon  modulus  requirement  Only  verdcal  bending  moment  composed  of  sdllwater 
and  wave  bending  moments,  is  considered.  The  inidal  yield  limit  state  is  expressed  as: 

g(X)  =  SM-Oy-Mj^-M^  (2.1) 

where  X  is  a  vector  of  the  random  variables,  ( SM,  Oy,  and  ),  and 
SM  is  the  secdon  modulus  amidship, 

Gy  is  the  yield  stress, 

is  the  sdllwater  bending  moment,and 
is  the  wave  bending  moment 

These  variables  are  taken  to  be  random  or  uncertain  and  are  assumed  to  be  stadsdcally 
independent. 
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12  Goieral  Characteristics  of  **  ABS  Ships" 


The  general  characteristics  of  several  ships  designed  to  the  minimum  requirements  of 
ABS  Rules  (including  minimum  section  modulus  requirements)  will  be  determined. 
These  ships  will  be  called  "ABS  Ships”.  Since  the  initial  yield  limit  state  is  the  only 
failure  mode  to  be  considered,  and  the  variables  in  Eq.  2.1  depend  only  on  L,  L/B,  and 
C^,  these  three  parameters  serve  as  the  factors  on  which  the  reliability  level  depends. 
They  are  specified  as  follows: 

L  :  from 91.Sm( 300 ft) to 366m (1200 ft) 

ITB  :  from  5.0  to  9.0 
C{)  :  from  0.60  to  0.85 

These  ranges  cover  most  ships  to  which  ABS  Rules  are  meant  to  apply.  The  value 
without  'tilde'  indicate  deterministic  characteristic  values. 


23  Strength  Considerations  of  "ABS  Ships" 

Because  of  variability  of  properties  of  steel  and  other  materials  used  in  marine 
structures  and  because  of  variability  in  production  and  fabrication  of  their  components, 
the  strength  of  identical  ships  will  not,  in  general,  be  identical.  In  addition,  uncertainties 
associated  with  residual  stresses  arising  from  welding,  the  presence  of  small  holes,  etc. 
may  affect  the  strength  of  the  ship.  These  limitations  and  uncertainties  indicate  that  a 
certain  variability  in  strength  or  hull  capacity  about  some  mean  value  will  result 

Additional  uncertainties  in  the  strength  will  arise  due  to  uncertainties  associated  with 
the  assumptions  and  methods  of  analysis  used  to  calculate  the  strength.  Further 
uncertainties  are  associated  with  possible  numerical  errors  in  the  analysis.  These  errors 
may  accumulate  in  one  direction  or  possibly  tend  to  cancel  each  other.  Whatever  the 
case,  the  above  uncertainties  have  to  be  reflected  in  any  reliability  or  failure  analysis. 


2.3.1  Section  Modulus 


Section  6  (Longitudinal  Strength)  of  ABS  Report  on  "Proposed  Change  to  Rules  for 
Building  and  Classing  Steel  Vessels"  September,  1991[1]  gives  the  minimum  required 
.ion  modulus  as  a  function  of  length  (L),  beam  (B),  and  block  coefficient  (C^  of  a 
ship  as  follows: 

SM  =  Ci-C2*L2-B*(Cb  +  0.7)  m-cm^ 
where  Cj  is  a  function  of  L,  and  C2  is  a  constant 


As  shown  in  Fig.  2.1,  the  section  modulus  is  assumed  to  be  lognormally  distributed 
with  a  coefficient  of  variation  of  4  %,  see  Ref.  6.  The  section  modulus  calculated  from 
the  ABS  rules  is  taken  as  the  mean  value. 


Lognormal  probability  density  function  (p.d.f.) 


Figure  2.1  Distribution  of  the  Section  Modulus. 


2.3.2  Yield  Strength 

The  yield  strength  distribution,  shown  in  Fig.  2.2,  is  assumed  to  be  lognormal  with  a 
coefficient  of  variation  of  7  %(Ref.  6),  and  with  a  mean  value  of  235  MPa  (34  ksi).  This 
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distribution  gives  a  probability  of  exceeding  ABS  permissible  stress  (17SMPa)  equal  to 
99.999%.  The  material  used  is  normal  strength  steel. 


Lognonnal  probability  density  function  (p.d.f.) 


Yield  Strength 


Figure  2.2  Distribution  of  the  Yield  Strength 


2.4  Loads  Applied  to  *'ABS  Ships" 

The  Stillwater  bending  monrent  was  obtained  from  the  1990  Rules[2],  the  latest 
available  at  the  time  the  work  was  conducted: 

Stillwater  Bending  Moment: 

Msw  =  10'^-Cst-L^-^'B-(Cb  +  0.5  )  kN-m  ( ‘90 ) 

Wave  Bending  Moment  Amidship  ( Sagging  Moment ): 

=  -ki*Ci‘L^*B*(  Cb  +  0.7  )  ‘lO’^  kN-m  (  proposed  for  '91 ) 

where  Cst,  k|,  are  constant,  and  is  a  function  of  L.  Hogging  moment  is  smaller,  and 
so  not  considered. 

Both  Stillwater  and  wave  moments  depend  on  length  (L),  beam  (B),  and  block 
coefficient  (Cb).  Fig.  2.3  shows  the  Stillwater,  wave,  and  total  bending  moment  variation 
with  ship  length  for  a  specified  block  coefficient  and  length-beam  ratio  as  an  example. 
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Appendix  1  shows  the  values  of  the  Stillwater  moment,  the  wave  moment,  the  ratio  of  the 
wave  to  Stillwater  moments  and  the  minimum  section  modulus,  all  calculated  according 
to  ABS  Rules  as  described  earlier  for  the  selected  ranges  of  length,  length  to  beam  ratio, 
and  block  coefficient. 

2.4.1  Stillwater  Bending  Moment  Distribution 

According  to  Soares  and  Moan[3],  tlw  Stillwater  bending  moment  fits  to  a  normal 
distribution.  In  this  investigation  it  is  assumed  that  the  value  given  by  ABS  is  the 
maximum  value  with  a  probability  of  exceedance  of  5  %.  The  large  variability  in  the 
Stillwater  bending  moment  calls  for  a  coefficient  of  variation  of  40%[3]  which  gives  the 
mean  value  of  the  distribution  to  be: 

1*  sw  ~  *^sw,ABS  (2-2) 

where  is  the  Stillwater  bending  moment  given  in  ABS  Rules  .  The 

distribution  is  shown  in  Fig.  2.4. 

Normal  Probablity  Density  Function  (p.d.f.) 


Figure  2.4.  Distribution  of  the  Still  Water  Bending  Moment 
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2.4.2  Wave  Bending  Moment  Distribution 


If  the  wave  loads  acting  on  a  mariiw  structure  can  be  represented  as  a  stationary 
Gaussian  process  (short-term  analysis),  then  at  least  four  methods  are  available  to  predict 
the  distribution  of  the  maximum  load.  These  methods  are  developed  for  application  to 
marine  structures  and  are  given  in  more  detail  in  [4].  In  this  report,  extreme  value 
distribution  based  on  upcrossing  analysis  [6]  is  used. 

The  wave  induced  bending  moment  given  by  ABS  is  modeled  as  an  extreme  value 
following  the  distribution  function[41: 


w 


Fw  ®*P  ) 


O 

w  ^  \y  21n  N 


(2.3) 


where  is  the  mean  of  the  distribution  and  is  the  standard  deviation.  N  is  the 
number  of  wave  bending  moment  peaks  and  is  the  mean  square  of  the  wave  bending 
moment  process.  The  value  given  by  ABS  is  assumed  to  be  the  mean  value  of  the 
distribution  [6],  and  Table  2.1  shows  how  the  coefficient  of  variation  varies  with  N. 
Choosing  N  to  be  1000,  which  is  equivalent  to  a  3  hour  storm  gives  a  coefficient  of 
variation  of  9  %.  Fig.  2.5  shows  the  distribution. 


N 

C.O.V. 

500 

10% 

1000 

9% 

2000 

8% 

Table  2.1 
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Extreme  value  probabiUly  density  function  (p.d.f.) 


Figure  2.5  Distribution  of  the  Extreme  Wave  Bending  Moment 

Appendix  2  gives  the  calculated  means  and  standard  deviations  of  the  Stillwater 
moment,  wave  moment,  and  the  section  modulus  according  to  the  distributions  described 
above  for  the  selected  ranges  of  L,  L7B  and  Cb- 

2.4.3  Comments  on  the  Ratio  of  Wave  to  Stillwater  Bending  Moments 
Given  by  ABS  Rules 

Inspection  of  the  calculated  values  of  Mg^,  M^,  and  M^/Mj^  according  to  ABS 
Rules  (Appendix  1),  leads  to  the  following  conclusions: 

1.  M  ^/Mg^  ratio  does  not  depend  on  L/B.  Hence,  M  ^/Mg^  can  be  written  as  a 
function  of  L  and  Cb  only. 

2.  Fig.  2.6  shows  the  ratio  M^/Mg^  as  a  function  of  L  for  two  extreme  values  of  Cb  (0.6 
and  0.85).  The  resulting  curves  are  more  or  less  parallel,  and  each  has  a  maximum  at 
L=152.5  m  and  a  minimum  at  L=366.0  m. 

3.  When  L  is  held  constant,  M  ^v^^sw  decreases  monotonically  as  Cb  increases. 


4.  As  a  result  of  the  above  observations,  all  Mg^/M^  values  fall  in  the  area  bounded  by 


Fig.  2.6  AfwlMsw  (  CbH>.6  Cb=0.85  ) 

as  a  function  of  length 
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2.5  Safety  Indices  and  Target  Reliability 


2.3.1  Reliability  Analysis  -  First  and  Second  Order 

The  reliability  analyses  are  carried  out  using  the  computer  program  CALREL  [3]  and 
fu’st  and  second  order  methods.  For  a  general  reference  of  these  methods  see  [6].  In  the 
reliability  analyses,  failure  is  defined  when  the  limit  sute  function,  g(20.  is  negative  ot 
zero.  X  is  it  vector  of  the  basic  random  variables.  i.e.  load,  material  and  geometrical 
properties.  After  transforming  the  basic  variables  into  standard  nomul  variates.^,  the 
program  determines  the  most  probable  failure  condition,  the  design  point,  through  an 
iterative  procedure.  The  design  point  has  the  cotHdinates  12*  where 

li*  =  -Pfl  (2.4) 


|3  is  the  safety  index  and  a  is  the  unit  row  vector  normal  to  the  tangent  plane  and  directed 
towards  the  failure  set,  see  Fig.  2.7.  FORM ,  the  First  Order  Reliability  Method,  replaces 
the  limit  state  surface,  g(X)  =  0,  with  a  tangent  hyperplane  at  the  design  point  in  the 
standard  normal  space  ,  while  SORM,  the  Second  Order  Reliability  Method,  replaces  the 
limit  state  surface  with  a  hypeq)araboloid  fitted  at  the  design  point  in  the  standard  normal 


Figure  2.7  The  First  and  Second  Order  Reliability  Methods 
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The  first  order  probability  of  failure.  Pf,  is  determined  from 

Pf»<I>(-P)  (2.5) 

where  <I>  is  the  standard  normal  distribution  function.  Fig.  2.8  shows  the  relation 
between  ^  and  Pf.  'P'  is  so  called  safety  or  reliability  index.  Tlw  higher  the  ^  value,  tlw 
lower  the  probability  of  failure,  and  the  higher  the  safety  margin  between  strength  and 
load.  The  reladonship  between  ^  and  Pf  given  in  Eq.  2.5  can  be  determined  numerkally 
from  the  properties  of  the  standard  normal  distribution  function  [15]. 

CALREL  was  used  to  calculate  reliability  indices  for  the  "ABS  ships"  covering  the 
entire  range  of  L,  L/B  and  described  earlier.  For  this  purpose,  the  limit  state  equation 
(2.1)  and  the  probability  distributions  given  in  sections  2.3.1, 2.3.2,  2.4.1,  and  2.4.2  were 
used  in  the  analysis.  Based  on  these  results  the  following  conclusions  are  made: 

1.  Holding  L,  L/B  fixed,  and  varying  €5  from  0.6  to  0.85 

As  shown  in  Fig  2.9,  the  safety  index  (^)  decreases  monotonically  as  the  block 
coefficient  increases. 

2.  Holding  L,  C5  fixed,  and  varying  L/B  from  5.0  to  9.0 

Fig  2.10  shows  that  ^  is  almost  constant  It  suggests  that  the  impact  of  L/B  on  ^  can 
be  neglected. 

3.  Range  of  P  for  different  L 

From  observations  1  and  2  above,  we  can  conclude  that  within  our  dimensions,  P 
varies  between  the  two  parallel  lines  shown  in  Fig.  2.11, which  shows  the  relation 
between  P  and  L  for  the  two  extreme  cases  (C^  =  0.6  and  0.85).  It  is  also  seen 
that  these  lines  have  the  same  pattern  as  lines  in  Fig.2.6.  Fig.  2.12  and  Fig. 

2.13  are  plotted  to  illustrate  the  relation  between  P  and  The  two  lines 

representing  the  boundaries  of  the  safety  indices  in  Figs.  2.12  and  2.13  are  plotted 
again  in  Fig.  2.14,  which  shows  that  they  fall  on  each  other.  This  suggests  that  p  can 
be  treated  as  a  function  of  Mw/Msw  only. 

4.  Table  2.2  shows  the  upper  and  lower  bounds  of  P  for  ship  length  varying  from 
152.5m  to  366m.  P  ranges  from  3.0236  to  3.3276  (see  also  Fig.  2.14),  and  its  average 
is  3.1918. 
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Urn) 

n 

r 

9U 

0.60 

3.2434 

32434 

0.85 

3.1635 

3.1635 

122.0 

0.60 

3.2953 

33070 

0.85 

3.2165 

32165 

152.5 

0.60 

3J276 

3.3272 

0.85 

32490 

32489 

183.0 

0.60 

33200 

33200 

0.85 

32416 

32416 

213.5 

0.60 

32933 

32933 

0.85 

32143 

32143 

244.0 

0.60 

32148 

32147 

0.85 

3.1343 

3.1343 

274.5 

0.60 

3.1992 

3.1992 

0.85 

3.1185 

3.1185 

305.5 

0.60 

3.1774 

3.1774 

0.85 

3.0962 

3.0962 

355.5 

0.60 

3.1389 

3.1389 

0.85 

3.0571 

3.0571 

366.0 

0.60 

3.1060 

3.1060 

0.85 

3.0236 

3.0236 

Table  2.2  Safety  Indices  of  ABS  Ships 
The  safety  check  equation  used  in  the  calculations  of  P  is  given  by  Eq.  2. 1 
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2.6  Comments  on  ABS  Rules  Regarding  Ship  Section  Modulus  Calculation 

The  following  conclusions  can  be  drawn  based  on  the  results  obtained  in  section 
2.5.1: 

1.  Safety  implied  in  ABS  Rules  for  longitudinal  strength  is  very  consistent  because  P 
varies  within  a  very  small  range.  However,  the  corresponding  ratio  of  the  upper  and 
lower  values  of  probability  of  failure  is  2.85.  This  means  that  some  room  for 
improvement  still  exists. 

2.  The  safety  index  depends  only  on  the  ratio  of  wave  bending  moment  to  Stillwater 

bending  moment.  This  makes  the  calibration  procedure  easier. 

3.  The  target  reliability  level  is  set  to  be  P  =  3.20,  which  is  approximately  the  average 
value  of  P  determined  earlier  for  the  ”ABS  Ships". 
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Probability  of  Fbilure  vs.  Safety  Index, 


^  'xapni 


Fig.2.8  Probability  of  Failure  versus  Safety  Index 
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lor  10  lor  lor*  lor* 

Probability  of  Fbilure 


I  I  I  I  I  j  T  I  \  I  — j  I  I - 1  I  I  T—  1  » - - 1  1 - T  i - 1  I 

0.5  0.6  0.7  0.8  0.9  1.0 

Block  Gxfncient  ( Q> ) 

Fig,  2,9  Safety  Index  versus  Cb  (L^l22mJUB^,0) 
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5.0  6.0  7.0  8.0  9.0  10.0 


Length'Beam  Ratio  ( L/B ) 

Fig,  2.10  Safety  Index  versus  UB  (  LsJ22m,Cb=0.6) 
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Fig,  2.11  Safety  Index  (  ) 

as  a  function  of  length 
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3.1 


1.5  1.54  1.58  1.62  1.66  1.7 

Mw/M$w . 

Fig.  2.12  Safty  index  versus  Mw/Msw  {LIB^S.0,Cb^.6) 


22 


Fig.  2.13  Safety  index  venus  Mw/Msw  (  UBssS.0,Cb^.85  ) 
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Safety  Index 


Fig,  2,14  Softy  index  versus  MwIMsw  (UBsS.O) 
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3.0  Calibration  Procedure 

Safety  factors  such  as  those  applied  to  yield  strength  and  to  loads  are  an  essential  part 
of  the  design  process.  In  the  probabilistic  methods,  this  need  resulted  in  the  introduction 
of  partial  safety  factors.  The  cumulative  effect  of  those  factors  is  such  that  the  resulting 
design  will  have  a  certain  reliability  level.  Thus,  code  developers  and  classification 
societies  may  determine  these  partial  safety  factors  that  ensure  that  the  resulting  design 
will  have  a  specified  reliability  level.  The  method  of  determining  these  partial  safety 
factors  for  a  given  safety  index  is  discussed  in  Reference[6]. 

The  objective  of  this  section  is  to  determine  partial  safety  factors  such  that  when 
applied  to  the  characteristic  values  of  Stillwater  moment,  the  wave  moment  and  yield 
strength,  the  resulting  hull  girder  section  moduli  for  all  ship  sizes  produce  constant 
reliability  index  equal  to  the  target  reliability  determined  earlier,  i.e.,  ^targen^-^- 
value  is  an  average  value  of  the  computered  safety  indices  for  the  ABS  ships  and  is 
selected  as  target  reliability  for  illustrative  purposes  only. 

3.1  Procedure  of  Calculating  Partial  Safety  Factors  for  "ABS  Ships'* 


As  described  above,  partial  safety  factors  are  used  in  the  calibration  procedure  to 
assure  a  specified  reliability  level.  For  the  current  case. 


SM  = 


at 


(3.1) 


where  •  Yw>  ^y  ^  partial  safety  factors  for  the  characteristic  values 
M^,  Oy  respectively. 


The  following  procedure  is  used  to  determine  the  partial  safety  factors  for  the  "ABS 
Ships" : 


1.  By  trial  and  error  determine  Y’s  and  0  in  Eq.  3.1  that  gives  the  ^target- 

2.  Find  out  for  different  ratios  of  the  value  of  3  determined  from  FORM  (or 

SORM)  using  the  '^s  and  <|>  obtained  in  the  first  step,  and  check  if: 

a.  the  obtained  P's  are  close  to  the  target  P,  and 

b.  the  obtained  Prange  is  smaller  than  that  of  ABS  rules. 
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3.  If  the  determined  Ys  ^  ^  give  ^'s  close  to  ^target  grange  is  smaller,  then  they 
can  be  used  in  the  new  calibrated  code,  otherwise  make  changes  in  them  to  satisfy 
the  two  criteria  a.  and  b.  above. 


3,2  Redesign  of  "  ABS  Ships”  and  Resulting  Safety  Indices 


The  procedure  described  above  can  be  implemented  as  follows.  Eq.  3.1  can  be 
rewritten  as: 


S_M 

Msw  <t>y<Jy 


(3.2) 


where  m  is  the  ratio  of  wave  bending  moment  to  Stillwater  bending  moment. 


It  is  obvious  that  in  Eq.  3.2  4)y  is  arbitrary,  so  we  set  it  to  be  0.86,  i.e.  a  material  or 
strength  safety  factor  of  I.IS.  Therefore,  if  we  can  find  two  ships  with  safety  indices 
equal  to  3.20,  a  pair  of  tentative  values  for  and  Yw  ^  determined.  One  ship  can 
be  directly  chosen  from  Table  2.2;  it  is  the  ship  with  Ls274.5m,  €5^.6,  and  ps3.i992. 
By  trial  and  error,  another  ship  can  be  found  by  changing  section  modulus  of  the  ship 
with  L=213.5m,  Cb=0.85  from  166690m-cm2  to  166374m-cm2  to  make  ^  equal  to 
3.2001.  The  values  of  7$^  7w  obtained  by  solving  the  resulting  two  equations 
when  the  values  are  substituted  in  Eq.  3.2.  The  resulting  ys  are: 


Ysw  =  1.103 
Yw  =115. 

Using  these  partial  safety  factors,  we  can  calculate  new  set  of  section  moduli  for 
which  we  perform  reliability  analysis  (CALREL)  to  determine  the  safety  index  for  every 
ship.  The  result  is  listed  in  Table  3.1  and  is  also  plotted  in  Fig.  3.1.  The  P's  in  Fig.  3.1 
are  very  close  to  each  other  (3.1980  <  p  <  3.2022),  as  compared  to  the  range  of  P  derived 
from  ABS  Rules.  Therefore,  the  calibrated  model  for  the  section  modulus  that  gives 
uniform  safety  for  all  ship  sizes  is  given  by  Eq.  3.1  with 

*ysw  =  1.103 
Yw  =  1.15 
()>  =  0.86. 
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Urn) 

Cs 

KSEB9 

9U 

0.60 

3.1999 

0.85 

3.2012 

122.0 

0.60 

3.1988 

0.83 

3.2004 

1S2.5 

0.60 

3.1980 

0.83 

3.1998 

183.0 

0.60 

3.1982 

0.83 

3.2000 

213J 

0.60 

3.1989 

0.83 

3.2001 

244.0 

0.60 

32003 

0.85 

3.2015 

274.3 

0.60 

3.1992 

0.83 

32017 

305J 

0.60 

32010 

0.83 

32018 

333J 

0.60 

32013 

0.83 

32020 

366.0 

0.60 

32018 

0.85 

3.2022 

Table  3. 1  Safety  Indices  of  Redesigned  ABS  Ships 
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Fig.  3.1  Safety4ndex  versus  Ship  Length 
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3.3  Benefits  of  the  Calibration 

The  main  benefit  that  accrues  from  the  redesign  exercise  according  to  the  new  safety 
check  format  is  uniform  reliability  and  structural  safety  among  different  ship  sizes, 
whichiii  some  cases  could  lead  to  weight  savings.  Code  calibration  exercises  such  as  this 
can  highlight  sometimes  large  differences  in  implicit  safety  levels  for  different  failure 
modes  in  a  structure,  a  situation  that  can  be  rectified  in  a  new  generation  reliability  based 
code. 
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4.  Development  of  Limit  States  for  an  Example  Ship 


As  stated  earlier,  the  objective  of  this  part  of  the  study  is  to  demonstrate  how  to  use 
reliability  technology  to  assess  the  level  of  risk  associated  with  an  existing  ship  or  with  a 
"drawing  board"  design.  For  this  purpose  an  existing  tanker  was  selected  as  an  example 
in  consultation  with  the  Project  Technical  Committee. 

Several  limit  states  are  formulated  and  applied  to  the  example  ship.  These  are:  the 
ultimate  limit  states  (deck  yielding,  fully  plastic  collapse,  and  instability  collapse),  the 
serviceability  limit  state  (local  buckling),  and  the  fatigue  limit  state  for  one  point  in  the 
deck.  Because  the  maximum  Stillwater  bending  moment  of  the  example  ship  occurs  in 
sagging  condition,  only  this  condition  is  considered  for  the  ultimate  and  serviceability 
limit  states.  IDetails  of  all  calculations  are  given  in  Appendices  3  through  7. 

4.1  Selection  of  the  Example  Ship 

A  tanker  designed  according  to  ABS  Rules  is  selected  as  the  example  ship.  The  main 
characteristics  are: 


Displacement 

149,000  tonnes 

L.O.A 

273.0  m. 

(895.1  ft) 

L.B.P 

260.0  m 

( 852.5  ft ) 

Beam 

42.0  m 

(  137.7  ft) 

Depth 

23.5  m 

(  77.0  ft) 

Draft 

16.0  m 

(  52.5  ft ) 

Cb 

0.710 

C  ^ 

The  elastic  section  modulus  at  deck  is  4.657675*10  m-cm  (236,851  in^-ft).  The 
nominal  yield  strength  of  the  material  used  is  259  MPa  (37.4  ksi). 


4.2  Formulation  of  Limit  States 

As  mentioned  earlier  the  limit  states  considered  in  this  demonstration  are: 


1.  Ultimate  strength  limit  state 


2.  Serviceability  limit  state 

3.  Fatigue  limit  state 

For  ships,  ultimate  limit  states  can  be  decomposed  into  two  modes  of  failure: 

a.  Failure  due  to  spread  of  plastic  deformation,  as  can  be  predicted  by  plastic  limit 
analysis  and  fully  plastic  moment  ( initial  yield  and  shake  down  moments  can  be  also 
classified  under  this  category )  [6]. 

b.  Failure  due  to  instability  or  buckling  of  longitudinal  stiffeners  (  flexural  or  tripping  ) 
or  overall  buckling  of  transverse  and  longitudinal  stiffeners  of  grillage. 

Serviceability  limit  states  are  associated  with  constraints  on  the  ship  in  terms  of 
functional  requirements  such  as  maximum  deflection  of  a  member  or  critical  buckling 
loads  that  cause  elastic  buckling  of  a  plate. 

Fatigue  limit  states  are  associated  with  the  damaging  effect  of  repeated  loading  which 
may  lead  to  loss  of  a  specific  function  or  to  ultimate  collapse.  This  particular  limit  state 
requires  an  independent  type  of  analysis. 

4.2.1  Ultimate  Strength  Limit  States 

Three  failure  modes  due  to  the  combined  action  of  wave  and  Stillwater  bending 
moment  are  considered.  The  ultimate  limit  state  can  be  described  as: 

Mu-Ms^-M^<0  (4.1) 


where 

M^  is  the  ultimate  hull  girder  moment  capacity  as  determined  by  the  critical  stress  of  the 
respective  failure  mode  and  the  effective  section  modulus. 

is  the  still-water  bending  moment. 

Mw  is  the  wave  bending  moment. 

M^  is  determined  for  each  failure  mode  as  follows: 

Deck  Initial  Yield 

Because  buckling  of  the  plates  in  the  deck  occurs  before  the  deck  initial  yield,  the 
effective  section  modulus  after  buckling  is  applied.  The  ratio  of  the  effective  section 
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modulus  to  the  elastic  section  modulus  is  calculated  to  be  0.98  (see  3.3  of  Appendix  3). 
The  critical  stress  is  then  the  material  yield  strength: 

SMeff  =  4.57-10^  m-cm^ 

Oct  =  259  MPa 


FuMv  Plastic  Collapse 

The  plastic  section  modulus  for  the  example  ship  is  calculated  according  to  [7],  and 
the  critical  stress  is  the  material  yield  strength.  The  details  of  the  calculations  are  given 
in  3.1  of  Appendix  3. 

SMp  =  5.8376*10^  m-cm^ 

Oct  =  259  MPa 
=  <'y 

Buckling  Instability 

The  elastic  section  modulus  is  used  and  the  critical  stress  is  the  buckling  stress  found 
by  applying  the  approximate  equations  described  in  [8].  These  equations  are  based  on 
beam  and  plate  theories  for  elastic  and  plastic  buckling.  The  elastic  section  modulus  of 
the  tanker  at  deck  is: 

SMp  =  4.65767-10^  m-cm^ 

and  the  critical  stress  due  to  buckling  depends  on  the  buckling  mode  as  follows: 
a.  Plates  between  stiffeners 

The  plates  between  the  longitudinal  stiffeners  are  considered  as  simply  supported 
isotropic  plates  under  uniaxial  compressive  load.  The  plate  collapse  stress  is  (see  3.2  of 
Appendix  3): 

Oct  =  238  MPa  (  ^  =0.92  ) 

Oy 
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b.  Stiffeners  and  effective  plating 


For  column  buckling  of  longitudinal  stiffeners  only  the  ultimate  limit  state  is 
considered  because  when  a  column  buckles  it  reaches  its  ultimate  strength  immediately. 
The  effective  plating  is  determined  from  buckling  considerations  since  the  plate  is  under 
edge  compression.  The  calculations  shown  in  3.2  of  Appendix  3  give  a  critical  stress  for 
pure  flexural  buckling  as: 

Oct  =  248  MPa  =0.958) 

Oy 

However,  coupled  tonional/flexural  buckling  stress  must  be  also  checked.  For  the 
example  tanker,  deck  longitudinal  stiffeners  have  a  single  plane  of  symmetry  which 
means  that  the  ultimate  limit  state  is  probably  governed  by  a  combination  of  torsional 
and  flexural  buckling.  For  this  condition,  the  critical  stress  is  (see  3.2  of  Appendix  3): 

=  170  MPa  (  ^  =  0.656  ) 

Oy 

c.  Cross-stiffened  panels 

Buckling  of  an  entire  stiffened  panel,  including  both  longitudinal  and  transverse 
stiffeners  is  considered  assuming  uniaxial  compressive  load.  A  panel  between  transverse 
and  longitudinal  bulkheads  is  shown  in  section  3.2  of  Appendix  3  together  with  the 
buckling  stress  calculations  according  to  reference[8].  The  resulting  critical  buckling 
stress  for  the  entire  panel  is 

Ocr  =  259  MPa 

d.  Summary,  Buckling  Limit  State  Strength 

Plate  between  stiffeners  238  MPa 

Flexural  buckling  of  stiffeners  248  MPa 

Tripping  of  stiffeners  170  MPa 

Cross  stiffened  panels  259  MPa 
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These  are  local  modes  of  failure.  The  ultimate  hull  girder  coUiqtse  monmtt  is 
calculated  in  item  e.  below. 

e.  Hull  Girder  Instability  Collapse 

In  the  1991  ISSC  proceedings,  report  of  the  Committee  on  Applied  Design[9],  the 
following  expression  was  used  for  the  {^proximate  determination  of  a  hull  girder 
instability  collapse  moment  in  sagging  condition: 

Ma  =  (-0.172+1.548^cp-0.368^cp^)*SMcOy 
is  the  compressive  strength  factm  given  by: 

=  (0.960+0.765X2+0.176B2+0.131X2b2+1.046X4)-0.5 

where 

X.  is  the  column  slenderness  of  a  critical  panel,and 
B  is  the  plate  slenderness  ratio. 

Appendix  4  shows  the  calculations  of  the  factor  for  the  example  tanker  and  the 
resulting  ultimate  moment  "My".  These  values  are 

♦cp  =  0.79  and 

My  =  0.82  SM^’Gy 

4.2.2  Serviceability  Limit  States 

The  serviceability  limit  state  can  be  expressed  in  the  same  form  as  for  the  ultimate 
limit  state: 

^serv. '  ^sw '  ^w  ®  (4*2) 


where 

M^rv.  moment  capacity  as  determined  by  the  critical  buckling  stress  in 

a  serviceability  limit  state. 
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A/ 

is  the  Stillwater  bending  moment 
is  the  wave  bending  moment 

The  critical  buckling  stress  of  local  plates  between  stiffeners  is  calculated  for  the 
example  ship  in  3.2  of  Appendix  3.  The  elastic  section  modulus  is  applied.  These  values 
are; 

SMg  =  4.65767*10^  m*cm^ 

Oct  =  227  MPa  (  ^  =0.870  ) 

Oy 


4.2.3  Fatigue  LimitJSiaK 

The  fatigue  limit  state  is  associated  with  the  damaging  effect  of  repeated  loading. 
There  are  two  approaches  to  the  fatigue  problem,  the  Palmgren-Miner  approach  based  on 
S-N  curves,  that  will  be  used  here,  and  the  firacture  mechanics  approach. 

The  S-N  curves  are  obtained  by  experiments  and  give  the  number  of  stress  cycles  to 
failure.  Such  curves  are  of  the  form: 

N*AS'"=  C  (4.3) 

where 

N  is  the  number  of  cycles  to  failure 
AS  is  the  stress  range 
m  is  the  inverse  slope  of  the  S-N  curve 
C  is  determined  from  the  S-N  curve  by 

logC  =  log  a  -2an,gN 

where 

a  is  a  constant  referring  to  the  mean  S-N  curve 
OiogN  ts  the  standard  deviation  of  logN 

The  fatigue  life  calculation  is  determined  based  on  the  assumption  of  linear 
cumulative  damage  (Palmgren-Miner  rule).  Application  of  this  assumption  implies  that 
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the  long-term  distribution  of  stress  range  is  replaced  by  a  stress  histogram  consisting  of 
an  equivalent  set  of  constant  amplitude  stress  range  blocks. 


The  time  to  failure  of  a  detail  can  be  expressed  as  [101  '• 


/V 

T» 


(4.5) 


where 

is  the  value  of  the  Palmgren-Miner  damage  index  at  failure. 

C  and  m  are  obtained  from  the  S-N  curves. 

B  is  the  ratio  between  actual  and  estimated  stress  range. 

is  a  stress  parameter. 


T,  Ap,  C  and  B  are  random  variables.  If  the  long-term  distribution  of  the  wave  process  is 
assumed  to  be  a  series  of  short-term  sea  states  that  are  stationary,  zero-mean,  Gaussian 
and  narrow  banded,  and  if,  in  addition,  the  structure  is  linear,  the  stress  range  will  follow 
a  Rayleigh  distribution  and  Q  is  determined  from[  10,11]: 


(m-l)/2  \n 

2n  ^  ‘21  Pj  ^oj 

j 


(4.6) 


where 

Pj  is  the  probability  of  occurrence  of  the  j-th  sea  state. 


Xpj ,  X2j  are  the  zero  and  second  stress  spectrum  moments  in  the  j-th  sea  state. 


respectively.  Note  that 
j-th  seastate. 


J-A  S  is 

2*  \/x<, 


the  frequency  of  the  stress  process  in  the 


The  fatigue  limit  state  function  is  expressed  as : 


(4.7) 


where  t  is  the  service  life  of  the  ship. 


5.  Developmeiit  of  Load  Modds  for  the  Example  Ship 

From  the  information  given  on  the  Tanker  example,  the  maximum  Stillwater  bending 
moment  is  1.9728*10^  kNm  and  it  occurs  in  sagging  condition.  The  maximum 
allowable  by  ABS  for  this  ship  is  3.022*10^  kNm. 

5.1  Wave  Bending  Moment  for  Ultimate  Limit  State 

The  r.m.s.  value  of  the  wave  induced  bending  moment  on  a  ship  can  be  estimated 
from  the  seakeeping  tables  in  [12].  Using  the  interpolation  procedure  described  in  that 
paper,  the  rms  of  the  bending  moment  can  be  determined  when  the  Froude  number,  the 
significant  wave  height  ,"H3",  the  beam/draft  ratio,  the  length/beam  ratio,  and  the  block 
coefficient  are  given.  Knowing  B/T,  L/B,  and  Cg  for  the  example  ship  and  assuming  the 
ship’s  speed  to  be 


12  knots  for  ^  3m 

8  knots  for  3m  <  ^  6m 

5  knots  for  6m  <  Hg. 

The  rms  of  the  wave  bending  moment  can  be  approximately  determined  for  any  sea  state. 

The  Wave  Bending  Moment  for  the  Ultimate  Limit  State 

For  the  ultimate  limit  state,  an  extreme  sea  condition  is  of  interest.  The  most  probable 
extreme  sea  condition  the  ship  is  likely  to  encounter  during  its  life  time  is  determined 
from  the  wave  data  along  its  route.  The  ship  is  assumed  to  remain  in  this  peak  sea 
condition  for  three  hours  (which  corresponds  to  N=1000  wave  peaks).  A  detailed 
procedure  for  this  short-term  analysis  is  described  in  reference[6].  The  wave  loads  in 
this  extreme  sea  condition  are  then  determined  and  the  corresponding  safety  indices  for 
the  ultimate  failure  modes  are  evaluated. 

Following  this  procedure  for  the  example  tanker,  the  rms  of  the  wave  bending 
moment  is  determined  for  a  significant  wave  height  of  12.2  m  (40  ft).  Section  5.1  of 
Appendix  5  shows  the  calculation  procedure.  The  resulting  rms  value  of  the  wave 
bending  moment  is 

=  rms  =  1.25398-10^  kNm  (5.1) 
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Assuming  that  the  wave  bending  moment  follows  the  same  distribution  as  described  in 
Section  2.4.2  with  N=10(X)  peaks,  the  mean  value  is  determined  by  Eq.  2.3  to  be 
4.855*10^  kNm.  For  comparison,  the  wave  bending  moment  given  by  1991  ABS  for  the 
example  ship  is  4.62*10^  kNm. 

Note  that  the  above  calculations  are  for  a  seastate  of  12.2  m  (40  ft)  wave  height  This 
particular  seastate  is  used  for  illustrative  purposes.  For  design,  a  storm  condition  with 
specified  return  period  should  be  selected  including  several  pairs  of  representative 
significant  wave  heights  and  characteristic  periods.  The  most  critical  ship  response  can 
be  thus  determined. 

5,2  Stress  Ranges  and  Number  of  Cycles  for  Fatigue  Limit  State 

The  sea  scatter  diagram  given  in  the  ISSC  proceedings[9]  and  shown  in  section  6.2  of 
Appendix  6  is  applied.  The  rms  value  for  every  sea  state  is  determined  and  the 
calculations  and  the  results  are  included  in  section  5.2  of  Appendix  5.  The  scatter 
diagram  used  is  for  the  Osebery  area  of  the  North  Sea. 
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6.  Reliability  and  Safety  Indices  of  the  Example  Ship 


In  this  section,  the  reliability  of  the  example  tanker  considering  both  the  ultimate  and 
fatigue  limit  states  is  determined.  Model  uncertainty  will  be  included  in  all  limit  state 
formulations  in  order  to  reflect  errors  resulting  from  assun^>dons  and  deficiencies  in 
analytical  or  en^irical  design  models  and  equations. 

6.1  Ultimate  Limit  States 


The  sagging  condition  is  considered  and  the  limit  state  is  expressed  as: 
g(X)  =  Xjj*SM*Cgy  -  ^w’^^s’^w 


(6.1) 


where 

SM  is  section  modulus. 

is  the  critical  failure  stress. 

is  the  sti'lwater  bending  moment 

is  the  wave  induced  bending  moment 

Xy  is  model  uncertainty  on  strength. 

is  uncertainty  in  the  model  of  predicting  the  Stillwater  bending  moment 

x^,  is  the  error  in  the  wave  bending  moment  due  to  linear  seakeeping  analysis. 

Xj  takes  into  account  nonlinearities  in  sagging. 

The  tilde  denotes  random  variables. 

The  distribution  of  model  uncertainty  parameters  are  shown  in  Table  6.1 


random  variable 

distribution 

mean 

c.o.v 

Xu 

N  (Normal) 

1.0 

0.15 

Xgw 

N 

1.0 

0.05 

Xw 

N 

0.9 

0.15 

Xc _ 

N 

1.15 

0.03 

Table  6.1  Distributions  of  Model  Uncertainty  Parameters 
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6.1.1  Deck  Initial  Yield 


Two  cases  of  the  Stillwater  bending  moment  are  consideied: 

In  CASE  1,  the  Stillwater  bending  moment  is  treated  as  a  deterministic  quantity  equal 
to  3.022-lO^cN  >m,  which  is  the  ABS  maximum  allowable  Stillwater  bending  moment 
for  this  ship.  The  effective  section  modulus  is  taken  as  the  meui  value.  Table  6.2  draws 
the  means  and  coefficients  of  variation  from  Itef.  [6]  of  the  random  variables  not  shown 
in  Table  6.1. 


random  variable 

distribution 

mean 

SM 

Lognormal 

4.57*10^  mcm^ 

0.04 

Orr 

Lognormal 

25.9  kN/cm^ 

0.07 

_ _ 

Extreme 

4.855*10®  kNm 

0.09 

Table  6.2  Distributions  of  Random  Variables  .CASE  1 

Appendix  7  shows  the  input/ouqiut  files  from  CALREL  printout.  The  safety  index  (P) 
equals  1.81,  which  implies  that  if  the  ship.while  loaded  at  its  iiuuumum  allowable  value 
of  the  Stillwater  bending  moment,  experiences  a  three  hour  storm  with  signiHcant  wave 
height  of  12.2m  (40  ft)  the  probability  of  failure  due  to  deck  yielding  is  Pf  »  3.5*10*^  for 
this  severe  storm. 

In  CASE  2,  the  Stillwater  bending  moment  is  treated  as  a  random  variable  with  mean 

equal  to  0.6*3.022*10^  to  be  consistent  with  Eq.  2.2.  Tables  6.1  and  6.3  give  the  random 

variables  and  their  distributions.  From  CALREL  for  this  case,  the  safety  index  (P)  equals 

.2 

2.25,  which  implies  a  probability  of  deck  yielding  of  Pf  »  1.2*10  . 

The  effect  of  correlation  between  the  Stillwater  bending  moment  and  the  wave 
bending  moment  is  investigated  next.  This  correlation  arises  because  of  a  weak 
dependence  of  the  wave  bending  moment  on  the  loading  condition.  CASE  2  is  repeated 
with  a  correlation  coefficient  of  0.2, 0.5,  and  0.8.  The  results  are  ^  2.23,  ^2.18,  and  ^ 
2.13,  respectively  for  this  severe  storm.  This  indicates  that  the  reliability  index  is  not 
very  sensitive  to  this  correlation  and  it  is  therefore  neglected  in  the  following  analyses. 
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random  variable 

distribution 

mean 

c.o.v 

SM 

Lognormal 

4.57*10^  mcm^ 

0.04 

S' 

Lognormal 

25.9  kN/cm^ 

0.07 

Normal 

1.813*10^  kNm 

0.40 

_ Mu. 

Extreme 

4.855*10®  kNm 

0.09 

Table  6.3.  Distributions  of  Random  Variables  .CASE  2 
6.1.2  FuUv  Plastic  Collapse 

The  random  variables  and  their  distributions  for  this  failure  mode  are  shown  in  Tables 
6.1  and  6.4.  The  limit  state  developed  in  Section  4.2.1  and  the  loads  determined  in 
Section  5  are  applied.  The  Stillwater  bending  moment  is  assumed  to  be  random.  This 
gives  a  reliability  ^3.15  and  a  probability  of  failure  of  8.3*10'^  for  the  severe  storm 
condition  considered. 


random  variable 

distribution 

mean 

c.o.v 

SM 

Lognormal 

5.838*10®  m-cm^ 

0.04 

Lognormal 

25.9  kNW 

0.07 

✓s/ 

Normal 

1.813*10®  kNm 

0.40 

Mu. 

Extreme 

4.855*10®  kNm 

0.09 

Table  6.4.  Distributions  of  Random  Variables,  Fully  Plastic  Colli4>se. 

6.1.3  Instability  Collapse 

Several  modes  of  failure  are  considered  under  instability  as  discussed  earlier.  These 
are: 

The  limit  state  developed  for  torsional/flexural  buckling  of  the  longitudinal  stiffeners 
is  applied  since  it  is  the  worst  mode  of  local  stability  failure.  The  load  is  as  determined  in 
Section  S,  and  the  Stillwater  bending  monwnt  is  assumed  random.  Tables  6.1  and  6.5 
give  the  random  variables  and  their  distributions.  From  CALREL,  P=0.S7  and  Pf  = 
2. 8*1  O' 1  for  the  severe  storm  condition  considered.  The  conditional  nature  of  this 


42 


probability  is  emphasized.  It  is  conditioned  on  encountering  this  severe  storm  condition, 
which  is  small.  The  mode  of  failure  is  also  local. 

The  hull  girder  instability  collapse  according  to  section  4.2.  l.d  is  considered  next 
This  gives  a  mean  value  of  =  212  MPa.  All  other  variables  remain  as  given  in  Table 
6.S.  The  resulting  safety  index  is  ^  s  1.49  and  Pf  ^  6.8*10*^,  again  conditional  on  the 
severe  storm  condition  considered. 


random  variable 

distribution 

mean 

SM 

Lognormal 

4.658*10^  m-cm^ 

0.04 

Orr 

Lognormal 

17.0  kNW 

0.07 

Normal 

1.813-10^  kNm 

0.40 

-■  ,  -  „ 

Extreme 

4.855*  10^*  kNm 

0.09 

Table  6.S.  Distributions  of  Random  Variables,  Instability  Collapse 


Fatigue  Limit  State 

Figure  6.1  shows  the  analyzed  detail,  which  is  a  welded  deck  longitudinal  to  the  deck. 
It  is  classified  as  class  D  according  to  classification  given  in  reference[13].  The  analysis 
is  concerned  with  one  fatigue  location.  No  system  aspects  are  considered.  The  limit 
state  function  is  given  as; 

x^  is  included  in  the  limit  state  as  a  modeling  uncertainty  to  take  into  account  the  error 
in  wave  bending  moment  prediction  using  linear  analysis.  The  other  variables  are  as 
described  in  Section  4.2.3.  The  stress  parameter,  calculated  in  section  6.1  of  Appendix  6, 
is  =  852  [  MN/m^]^[scc]'*  and  from  the  S-N  curve,  the  mean  value  of  C  =  1.52*10^^ 
MN/m^. 

The  analysis  is  performed  with  the  random  variables  distributed  as  shown  in  Table 
6.6.  The  reliability  index  P  equals  2.44,  and  the  probability  of  failure  is  7.3* lO'^  ovct  a 
lifetime  of  20  years. 
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longitudinal 

sulfcncr  dctt  plate 

(450  x  30)  1,21  mm 


deck  plate 


j. _ a 

4 

^lOOOmm^ 

- - 4 

^lOOOmm^ 
■  ■ 

L  450  X  30 

r 

A  -  A 


Figure  6.1  Detail  Considered  in  the  Fadgue  Analysis. 


random  variable 

distribudott 

mean 

Lognormal 

1.44 

0.1S 

6 

Lognormal 

1.5M0‘* 

0.40 

A/ 

B 

Lognormal 

1.02 

0.10 

*w 

Ntmnal 

0.90 

0.1S 

Table  6.6.  Distributions  of  Random  Variables,  Fadgue 
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6.3  Summary  of  Safety  Indices 


The  following  is  a  summary  of  the  calculated  probabilities  of  failure: 

a)  Deck  initial  yield  0.012  (Global) 

b)  Fully  plastic  condition  0.00083  (Global) 

c)  Instability  Stripping)  0.28  (Local) 

d)  Hull  girder  ultimate  moment  0.068  (Global) 

e)  Fatigue,  20  years  0.007  (Local) 

It  is  to  be  emphasized  that  these  values  are  conditional  on  the  severe  seastate  assumed, 
in  the  case  of  items  a)  through  d).  The  unconditional  probabilities  of  failure  are  expected 
to  be  lower  since  the  shown  values  in  items  “c"  and  ”d"  must  be  multiplied  by  the 
probability  of  encountering  the  severe  storm  condition  used  in  their  calculations.  The 
fatigue  (item  e)  is  unconditional  value  calculated  for  one  detail  over  the  20  year  life  of 
the  ship. 
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7.  Ternunology  Associated  with  Structural  Rdiability 

The  aim  of  this  chapter  is  to  define  the  terminology  associated  with  the  structural 
reliability  of  ships  and  offshore  structures.  The  following  are  considered; 

•  Load  terminology 

•  Strength  terminology 

•  Stmctural  reliability  terminology 

The  terminology  defined  addresses  those  ^rms  associated  with  probability,  statistics  and 
reliability  as  used  in  engineering. 

7.1  Load  Terminology 

The  following  terms  are  primarily  used  with  loads,  although  some  of  the  terminology 
is  more  general,  and  related  to  statistics  and  random  processes. 

Deterministic  Process 

If  an  experiment  is  performed  many  times  under  identical  conditions  and  the  records 
obtained  are  always  alike,  the  process  is  said  to  be  deterministic.  For  example,  sinusoidal 
or  predominantly  sinusoidal  time  history  of  a  measured  quantity  are  records  of  a 
deterministic  process. 

Random  Process 

If  the  experiment  is  performed  many  times  when  all  conditions  under  the  control  of 
the  experimenter  are  kept  the  same,  but  the  records  (usually  a  time  history)  continually 
differ  fK»n  one  another,  the  pocess  is  said  to  be  random.  The  degree  randmnness 
depends  on  (1)  understanding  of  the  factors  involved  in  the  experiment  results,  (2)  the 
altility  to  control  them.  The  outcome  ctf  a  raiKlom  process  at  any  given  instant  of  time  is 
a  random  variable.  Time  history  of  wave  elevation  and  strain  gage  records  taken  aboard 
a  ship  may  be  ccmsidered  as  random  {nrocesses. 

Random  Variable 

Different  values  of  a  random  variable  have  (Efferent  chances  (frequencies)  of 
occurrence.  A  random  variable  thus  has  a  probability  density  function.  Examples  of 
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random  variaUes  are  die  wave  ht»iMKng  moment,  die  still  water  bending  moment,  and 
material  yield  strength. 

Prohabtlitv  Denritv  Funetka 

The  probability  density  fiiiicdon  defines  the  relative  frequencies  of  occurrence  of  a 
random  vaxiaUe  (e.g.,  wave  height  or  wave  bending  moment).  The  function,  usually 
denoted  f(x),  where  X  is  die  random  variaUe,  has  the  following  properties: 


1)  The  probability  of  occurrence  of  fraction  of  the  random  variable  X  which  lies 
between  x  and  x-Kdx  is  f(x)dx,  Le., 

P[xSX^x  +  dx]  =  f(x)dx 

2)  The  probability  that  a  sample  of  the  variable  lies  between  a  and  b  is: 

P[a^XSb]  =  Jf(x)dx 

3)  The  probability  that  X  lies  between  -«-» and  +» is  unity. 

4)  P[x  =  a]  =  0  where  a  is  a  constant 


ProbabUitv  Distribution  Function 

Also  called  the  cumulative  distribution  function,  and  denoted  F(x),  this  defines  the 
probability  that  the  random  variable  X  is  less  than  or  equal  to  a  given  value  x,  i.e.. 
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F(x)*-|^f(x)dx 


EMccdcncc  Probability 

This  is  the  probability  that  a  random  variable  X  (e.g.,  wave  bending  moment) 
exceeds  a  specified  value  x,  and  is  given  in  terms  of  the  probability  distribution  function 
as  1  -  F(x),  since 


1-F(x)=£f(x)dx 
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Peroenrile 

PcKcemile  vi^nes  of  «  xandocB  ^raiiiMe  X  tre  those  ^r^aes  cocseqpoadii^  to  qKctfied 
values  of  the  cumulative  distribution  function  F(x).  A  50-pacentile  value  thus 
COTiesponds  to  x  such  that  F(x)  =  0.5.  This  particular  percentile  is  also  the  median  value 
of  the  random  variable.  A  9S-percentile  value  is  a  value  such  that  F(x)  =  0.9S,  i.e.,  only 
5%  of  the  outcomes  of  the  random  variable  are  expected  to  lie  above  it 


Mean.  Median  and  Mode 

For  a  given  probability  density  function  f(x)  relating  to  a  random  variable  X,  the 
mean  or  average  value  p.  is  given  by 


p  =  E(x)  =  j~  xf(x)dx 

where  E(x)  denotes  the  "expected  value"  of  X. 

The  median  value  of  X,  denoted  x,  is  defined  from  the  cumulative  distribution 
function  F(x)  as 


x  =  F'‘(0.5) 

i.e.,  it  is  a  value  of  X  corresponding  to  a  cumulative  distribution  function  of  0.5.  This 
implies  that,  on  the  average,  1/2  the  outcomes  of  the  random  variable  will  lie  below  x 
and  1/2  above  it. 
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The  mode  of  a  nmdom  vaiiaUe  X  is  the  value  of  X  cone^xmding  to  the  peak  the 
probability  density  for  the  random  variable.  The  mode  is  also  called  the  most  probable 
value  of  the  random  variable  (e.g.,  most  probable  wave  bending  moment). 


Mean  Square  Value 

The  mean  square  value  of  a  random  variable  X  is  defined  by 

E(x*)  =  £x*  f(x)dx 

and  its  root-mean-square  or  r.in.s.  value  is  simply 


Variance  and  Standard  Deviation 

The  variance  of  the  random  variable  X  is  defined  by 

=  E(x-pJ*  =£J(x-p,)*  f(x)dx  =E(x*)-p^ 

The  standard  deviation  of  the  random  variable  is  a.  The  standard  deviation  is  a  measure 
of  spread  of  the  random  variable  about  the  mean  value.  Note  that  for  a  zero  mean 
variable,  the  variance  and  the  mean  square  value  are  numerically  the  same.  This  is 
approximately  true  for  both  waves  and  wave  bending  moment  assuming  linear  first  order 
theory  holds. 
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tjoemcientot  VMianon 

The  coeCGdem  (tf  varutioa  8  a  nuKk)m  vtriabte  X  is  defined  by 

6=2^ 

11. 

where  a  and  ^  are  the  standard  deviation  and  the  mean  value.  The  coe£Bcient  ai 
variatitm  is  a  non-dimensional  measure  of  the  ^nead  of  the  random  variable  outcomes 
about  the  mean  value.  The  coefficimt  variation  of  wave  heights  and  wave  bending 
moments  over  a  long  period  of  rime  is  expected  to  be  hi^  (80-100%).  The  coefficient  of 
variation  of  the  extreme  values  of  diese  quantities  over  a  short  period  of  time  in  a  severe 
sea  state  is  much  smaller  (7-20%). 


Joint  ProbabUitv  Densitv  Function 

The  joint  probability  density  functicm  of  two  random  variables  x,  and  Xj  defines  the 
fiequency  of  mutual  occunroice  of  r«ro  random  variables  and  has  the  following 
properties: 

1)  P(x,  <X,  ^Xj*fdx,  fix,  <X2  ^Xj+dx,  ~f(x,,X2)dx,  dx2 

2)  P[a,  <  X,  ^  b,  n  a2  <  X2  ^  b2] = )  dx,  dx2 

3)  P[-«»<X,  <+«n-«»<X2  <+«»}=  f(x,,X2)dx,  dx2  =  1 

where  n  indicates  the  mutual  occurrence  (intersectitm)  of  two  events. 

A  related  joint  distribution  function  defining  cumulative  probabilities  may  also  be 
defined.  The  definitions  may  be  extended  to  nxne  than  two  random  variables. 

The  joint  density  and  distribution  functicms  for  random  variables  contain  the 
occurrence  probability  and  also  correlation  infCTmation. 


rnyimance 

The  covariance  of  two  random  variables,  Xj  and  X2  is  defined  as 

(‘x.o,  -  e{[x,-E(x,)I»,-E(x,)1} 


- 

where  and  are  the  means  of  the  individual  random  variables,  and  f(X|,  is  their 
joint  density  functicm. 

Indepeiident  Random  Variables 

Two  random  variables  Xj  and  X2  are  independent  if  their  joint  density  function  is 
equal  to  the  product  of  dieir  individual  denrities 

f(X,,X2)  =  f(X|)f(X2) 

where  f(X),  X2)  is  the  joint  density  function  and  fCx^)  and  f(x2)  are  the  individual  (also 
marginal)  density  functions.  The  outcomes  of  independent  random  variables  occur 
without  any  reference  to  one  another.  Nmmally  in  reliability  analysis,  strength  and  load 
are  considered  indqtoident  random  variables. 


Dependent  Random  Variables 

Two  random  variables  X,  and  X2  are  dq)endent  if  their  joint  density  function  is  not 
the  product  of  the  marginal  densities.  The  outcome  of  any  one  of  the  random  variables  is 
dependent  on  the  outcome  of  the  other,  i.e.,  there  is  a  correlation  between  the  realization 
of  one  random  variable  and  realizations  of  the  other.  For  Xj  dependent  on  X2,  the 
following  is  true: 


f(x,/x2)  = 


^X,Xx(*l  ^*2) 

f(X2> 


where  f(x,/x2)  is  the  conditional  density,  fCxj)  is  a  marginal  density,  and  /x^)  is 

the  joint  density  evaluated  with  X]  given  X2. 
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Bounded  Random  Variables 

The  definitions  of  probability  density  and  distrilnitiim  functions  given  in  this  secticm 
assume  that  random  variable  outcomes  lie  in  the  interval  -«>  <  X  <  4«o.  Here,  the  bounds 
on  the  random  variable  are  -<»  and  -h*.  For  sonae  random  variables,  the  upper  and/or 
lower  bounds  may  be  different.  For  example,  material  yield  strength  is  always  a  positive 
quantity,  and  its  lower  bound  is  zero.  An  upper  bound  on  a  load  is  sometimes  used 
resulting  in  a  truncated  probal^ty  density  function. 


Conclanon  Q, 


The  correlation  coefficient  for  two  random  variables  Xj  and  X2  is  defined  by 


*t*t 


where  is  the  covariance  of  x,  and  X2,  and  the  o  are  the  standard  deviations.  The 

correlation  coefficient  always  lies  between  -1  and  +1.  If  the  cmielation  coefficient  is 
zero,  the  variable  outcomes  are  tmoortelated.  The  correlation  coefficient  IS  a  first  order 
measure  of  dependence  between  outcomes  of  two  random  variables.  A  zero  correlation 
is  a  weaker  condition  than  independence.  Non-conelated  random  variables  are  not 
necessarily  independent,  but  independent  random  variables  are  necessarily  unctnrelated. 
Positive  correlation  means  that,  in  general,  if  the  outcomes  of  one  random  variable 
increase,  the  outcomes  of  the  other  will  also  increase.  Negative  correlation  means  that 
the  outcomes  will  generally  be  in  opposite  directions. 

The  wave  bending  moment  is  weakly  correlated  to  the  stiUwater  bending  moment 
since  both  depend  on  the  weight  drstributkxr  along  ship  length. 


Conditional  ProbabUitv  and  Baves  Theorem 

A  conditional  probability  is  denoted  P[A/B]  when  A  is  one  event  and  B  is  another 
event  on  whose  outcome  A  depends  on.  An  example  of  a  conditional  probability  is  a 
probability  of  structural  failure  calculated  for  a  given  sea  state.  The  actual  lifetime 
probability  of  failure  will  be  different  if  all  the  sea  states  are  considered.  Bayes' 
Theorem  applies  to  conditional  events.  By  Bayes'  Theorem,  the  probability  that  event  A 
occurs  conditioned  on  the  probability  that  event  B  has  already  occurred  is  given  by 


PiAIB)  = 


PjAnB) 

PiB) 
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where  A  and  B  are  the  evatt  domains  and  AnB  is  tbar  intenectk)ii«  Le..  the  outcome 
space  that  contains  both  A  and  B  at  the  same  time  (mutual  occuneaceX 

StttiPiugy  Random  Preccaa 

A  nmdom  process  is  statumaiy  if  die  probability  densiQr  function  of  its  outcomes 
does  not  dqiend  on  time,  Le.,  the  same  piobal^ty  density  function  is  obtained  for  an 
ensemble  of  realizations  of  the  random  process  at  any  ghrm  time  as  at  any  other  time. 
This  also  means  that  statistics  that  are  dependent  on  the  probability  density  function,  e.g., 
mean  and  mean-square  value,  are  also  independent  of  time.  The  second  Ofxler  (joint) 
probability  density  function  of  the  outcomes  at  two  instants  of  time  depoids  on  die  time 
lag  between  them  and  not  <»  each  individually.  Time  history  of  waves  or  wave  bending 
moment  are  usually  considered  stationary  over  a  short  period  of  time  (up  to  3  hours). 

Ereodic  Hypothesis 

This  states  that  a  single  sanqile  function  is  quite  typical  of  all  other  sanqile  functions 
representing  realization  of  a  random  process.  Therdbre  we  can  estimate  the  various 
statistics  of  interest  by  averaging  over  time  using  the  one  realization  rather  than 
averaging  over  an  ensemble  of  realizations.  An  ergodic  random  process  is  necessarily 
stationary.  A  stationary  random  process  is  not  necessarily  ergodic. 

Exteme  Value 

The  extreme  value  of  a  random  process  is  the  largest  value  over  a  period  of  time. 
Each  realization  of  the  random  process  will  have  an  extreme  value.  Thus  there  is  also  a 
distribution  of  extreme  values,  i.e.,  the  extreme  value  is  a  random  variable  that  has  its 
own  special  distribution,  mean  value,  variance,  etc.  One  may  speak,  therefore,  of 
extreme  value  distribution  of  wave  heights  or  wave  bending  moments. 

Most  Probable  Extreme  Value 

This  is  the  value  of  the  random  variable  corresponding  to  the  peak  of  the  extreme 
value  density  function,  i.e.,  the  mode.  Thus,  the  most  probable  extreme  wave  bending 
moment  is  the  mode  value  of  the  extreoM  bending  moment  density  function,  i.e.,  the  value 
of  the  moment  at  the  peak  of  the  density  function. 

Asymptotic  Distributions  of  the  Extreme  Value 

The  extreme  value  distribution  for  a  random  process  with  defined  probability 
characteristics  for  the  outcome  (e.g.,  a  Gaussian  random  process)  is  a  function  of  time,  or 
equivalently,  the  number  of  peaks  within  the  time.  As  time  or  number  of  peaks  iiKsease, 


55 


the  distribution  of  the  extreme  value  shifts  to  the  right.  The  asymptotic  distribution 
corresponds  to  an  infinite  length  of  time  or  number  of  peaks.  The  asymptotic  form  of  the 
extreme  value  distribution  depends  largely  on  the  tail  behavior  of  the  "initial"  distribution 
of  outcomes  of  the  random  process.  Gumbel  showed  that  the  asymptotic  distribution 
takes  one  of  three  forms:  a  double  exponential  form,  an  exponential  form  and  an 
exponential  form  with  an  upper  bound. 

Order  Statistics 

The  distribution  of  the  largest  peak  (e.g.,  largest  wave  bending  moment)  in  a 
sequence  of  N  peaks  of  a  random  process  can  be  determined  using  order  statistics, 
assuming  that  the  peaks  are  independent  and  identically  distributed.  The  cumulative 
distribution  function  of  the  largest  peak  is  given  by 

Fz^(Z)  =  P[max 

=  [F,(z,e)]N 

where  F2(z,e)  is  the  initial  cumulative  distribution  function  of  the  peaks  and  e  is  the 
spectral  bandwidth  parameter.  The  corresponding  probability  density  function  is  given 
by  differentiating  the  cumulative  distribution  function: 

=  N[F,(z,e)]N->  •  f,(z,E) 
where  f2(z,E)  is  the  initial  p.d.f.  of  the  peaks. 

Expected  Maximum  Value: 

The  expected  value  (average)  of  the  maximum  peak  (e.g.,  wave  bending  moment)  in 
a  sequence  of  N  peaks  of  a  zero  mean  Gaussian  random  process  was  determined  by 
Cartwright  and  Longuet-Higgins,  and  is  approximated  by 


where  C  =  0.5772  =  Euler's  constant.  Here,  is  the  area  under  the  power  spectral 
density,  i.e.,  the  mean  square  value  of  the  process. 

It  should  be  noted  that  the  most  probable  extreme  value  (i.e.,  the  mode)  is  given  by 
the  above  equation,  but  with  the  second  term  on  the  right  hand  side  deleted. 
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Nanpw  Band  Process 

This  is  a  random  process  whose  time  lealizatitHis  are  such  that  there  is  one  peak 
between  every  upcrossing  and  every  downcrossing  of  the  mean  level.  Process  "cycles” 
are  thus  discernible.  The  power  q)ectral  density  functimi  of  the  process  realization  has  a 
central  tendency,  i.e.,  it  is  clustered  about  a  central  frequency.  The  peaks  of  a  zero  mean 
narrow  band  Gaussian  random  process  have  the  Rayleigh  distribution  function  given  by 

where  mg  is  the  mean  square  value  of  the  process,  also  equal  to  the  area  under  the  energy 
spectram  for  the  process. 

Records  of  waves  and  wave  bending  moments  over  a  shon  period  of  time  (3  hours) 
are  usually  considered  to  be  narrow-band  processes. 

Average  of  Highest  1/m-th  Value 

This  is  the  average  value  of  the  highest  1/m-th  peaks  in  a  random  process.  For  a 
random  process  whose  peaks  are  Rayleigh  distributed. 

Average  of  1/3  highest  values  *  2 

Average  of  1/10  highest  values  =  2.5S 

Average  of  1/1000  highest  values  =  3.85 

where  m^  is  the  mean  square  value  of  the  process.  The  multipliers  shown  are  for 
amplitudes  rather  than  heights  (double  amplitudes).  The  average  of  1/3  highest  values  is 
also  called  the  significant  value.  These  multipliers  may  be  used  for  waves  and  wave 
beading  mmnents  andmay  err  slightly  cm  the  ccmservative  side. 

7.2  Strength 

The  following  terms  related  to  strength  are  now  defined:  failure  mcxles,  limit  state 
function,  and  ultimate,  serviceability  and  fatigue  limit  states.  Limit  state  exceedence 
probability  is  then  defined,  and  contrasted  to  the  probability  of  failure.  Also  in  this 
section,  terminology  related  to  the  classification  of  uncertainties  is  given.  Some  of  this 
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tenniiiology  is  general,  but  their  use  is  relevant  to  strength  variability,  a 
with  strength  parameters.  System  failure  modeling  is  also  considered  in  this 

Failure  Mode: 

A  failure  mode  refers  to  a  particular  physical  mechanism  which  a  structure  or  a 
part  of  it  fails.  Failure  modes  for  ships  address  plastification.  buckling,  fatigue  and 
fracture.  As  an  example,  buckling  failure  modes  include  plate  budding,  stiffener  flexural 
buckling,  stiffener  tripping,  and  overall  buckling  of  the  gross  panel. 

Ultimate  Limit  State: 

The  ultimate  limit  state  considers  structural  performance  or  safety  mar^  under 
extreme  (typically  lifetime  maximum)  loads.  The  ultimate  limit  state  can  be  fiuther 
decomposed  into  two  modes  of  failure: 

a.  Failure  due  to  spread  of  plastic  deformation.  e.g.,  as  predicted  for  beams  by  plastic 
limit  analysis.  The  initim  yield  moment  for  a  beam  can  also  be  classified  under 
this  category. 

b.  Failure  due  to  instability  or  budding.  e.g..  of  panel  longitudinal  stiffeners  in  the 
flexural  and  tripping  modes,  or  the  overall  "grillage"  buckling  of  a  gross  panel 
consisting  of  longitudinal  and  transverse  stiffeners. 

SetviceabiHtv  Limit  State: 

The  serviceability  limit  states  are  associated  with  constraints  on  the  marine  structure 
in  terms  of  functional  requirements  such  as  the  maximum  deflection  of  a  member  or 
critical  buckling  loads  that  cause  elastic  buckling  of  plating. 

Fadeue  Limit  State: 

The  fadgue  limit  state  is  associated  with  the  damaging  effect  of  repeated  loading 
which  may  lead  to  a  loss  of  specific  funcdon  or  to  uldmate  collapse.  Fadgue  limit  state 
capacity  for  structural  details  is  typically  defined  using  S-N  ciuves.  while  the  demand  is 
defined  in  terms  of  the  lifetime  stress  range  versus  number  of  cycles  histogram. 

Limit  State  Funcdon: 

This  is  a  function,  often  denoted  G(2D  where  X  is  a  vector  of  basic  variables,  that 
characterizes  the  safety  margin  in  a  given  mode  of  failure.  A  simple  limit  state  function 
may  be 
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G(Oy,  O)  =  Oy  -  O 


where  Oy  is  the  yield  strength  of  the  material,  and  o  is  the  load  effect  (stress).  Note  that 
limit  state  exceedence  ("failure")  implies 


G  ^  0 

Limit  state  functions  are  traditionally  formulated  in  this  capacity  minus  demand  fmm. 
The  basic  variables  in  the  limit  state  equation  are  random  because  of  inherent  variability 
or  model  uncertainties. 

Limit  State  Exceedence  Probability 

The  probability  of  reaching  or  exceeding  a  specified  limit  state  is  determined  from 

F 

where  f^(20  is  the  joint  probability  density  function  of  the  basic  variable  vector  X  The 
domain  of  integration  F  is  over  the  unsafe  region  of  the  limit  state  function  where 
demand  exceeds  capability.  The  integral  is  multi-fold.  In  terms  of  the  limit  state 
equation,  the  domain  of  integration  is  defined  by  G(20  ^  0.  To  the  extent  a  limit  state 
equation  may  address  local  phenomena,  e.g.,  yield  at  a  point,  serviceability,  e.g., 
deflections,  etc.  in  addition  to  catastrophic  events,  interpreting  the  limit  state  exceedence 
probability  as  the  probability  of  "failiKc”  of  the  structure  should  be  done  with  care. 

It  should  also  be  noted  that  limit  state  exceedence  probabilities  calculated  are  c^ten 
conditional  on  certain  environmental  c  rents,  e.g.,  occurrence  of  a  certain  severe  storm. 

Probability  of  Failure 

Although  actuarially  speaking,  this  should  refer  to  the  probability  that  the  structure 
catastrophically  fails,  the  term  is  generally  and  widely  used  as  a  substitute  for  limit  state 
exceedence  probability,  i.e.,  the  probability  that  the  demand  exceeds  the  capability  in  any 
given  limit  state  (including  exceedence  of  deflection  and  clastic  buckling  stress). 

Uncenaintv  Classification 

Uncertainties  which  contribute  to  the  variability  of  physical  strength  parameters  may 
be  classified  as 


59 


•  inherent  uncertainties 

•  model  uncertainties 

They  may  also  be  classified  as  subjective  and  objective  uncertainties.  The  classifications 
while  illustrated  here  with  stren^  parameters,  are  also  relevant  to  loads  and  load 

_ j  -1  - 

moQCxS* 

Objective  Uncertainties 

These  are  uncertainties  associated  with  random  variables  for  which  statistical  data  can 
be  collected  and  examined.  They  can  be  quantified  by  a  mean,  a  coefficient  of  variation, 
and  a  form  of  the  probability  distributicm  function  derived  from  available  statistical 
information.  The  variability  in  the  yield  strength  of  steel  is  an  example. 

Subjective  Uncertainties 

These  are  uncertainties  associated  with  the  lack  of  information  and  knowledge.  They 
are  typically  quantified  on  the  basis  of  the  engineer's  prior  experience  and  judganent 
Exanq)les  of  these  include  assumptions  in  the  analysis,  error  in  the  design  model,  and 
empirical  formulae.  The  following  subjective  uncertainties  contribute  to  strength 
variability: 

a)  Effectiveness  of  plating,  e.g.,  due  to  shear  lag 

b)  Use  of  Navier  hypothesis  in  calculating  hull  girder  response 

c)  Initial  deformation  and  residual  stress  effects 

Inherent  Uncertainties 

This  kind  of  uncertainty  is  inherent  to  the  variable,  and  cannot  be  reduced  because  of 
additional  information.  This  is  a  term  timt  in  ^c^y  cases  may  involve  the  same  sources 
as  "objective"  uncertainties.  Examples  are  the  inherent  variability  of  wave  heights, 
extreme  w!<ve  bending  moment  or  the  variability  in  yield  strength. 

Model  Uncertainties 

These  uncertainties  arise  because  of  otots  in  the  prediction  models  as  they  represent 
reality.  They  can  be  reduced  with  additional  information.  Model  uncertainties  are 
typically  estimated  based  on  comparing  the  analysis  procedure  with  experimental  data,  or 
in  some  cases  using  professional  judgement  or  other  indirect  information  such  as  the  non¬ 
occurrence  of  cracks  in  relation  to  expectation.  Some  sources  of  model  uncertainties  are 
described  under  "subjective  uncertainties".  The  largest  model  uncertainty  in  marine 
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structures  usually  relates  to  loads  such  as  slamming  loads.  Strength  prediction 
techniques  (e.g.,  for  buckling  strength)  also  have  their  own  model  uncertainties.  This 
type  of  uncertainty  is  usually  quantified  in  terms  of  a  bias  (i.e.,  actual  value  to  predicted 
value  ratio)  and  a  coefficient  of  variation. 

Structural  System  Modeling 

The  behavior  of  a  structure  that  can  fail  in  mote  than  one  mode  of  failure  is  modeled 
for  structural  reliability  evaluation  purposes  using  structured  representations  of  system 
behavior.  Series,  parallel  or  general  system  representations  are  usual.  A  general  system 
representation  may  take  the  form  of  a  cut  set  (parallel  subsystems  connected  in  series) 
representation  or  a  link  set  (series  subsystems  connected  in  parallel)  representation. 
Failure  tree  representations  are  also  possible.  Reference  is  made  to  [6]. 

Series  System: 

A  series  system  is  one  that  is  composed  of  links  connected  in  series  such  that  failure 
of  any  one  or  more  of  these  links  constitute  a  failure  of  the  system,  i.e.,  "weakest  link" 
system.  In  the  case  of  the  primary  behavior  of  a  ship  hull,  for  example,  occurrence  of 
any  one  of  a  number  of  modes  of  failure  will  constitute  failure  of  the  hull.  The  multiple 
failure  modes  can  then  be  modeled  as  a  series  systerrL 

Parallel  System 

In  a  parallel  system,  all  links  along  the  failure  path  must  fail  for  the  structure  to  fail. 
An  example  is  a  multicomponent  redundant  structure  such  as  a  fixed  offshore  platform, 
in  which  a  failure  path  is  the  failure  of  a  group  of  members  which  leads  to  system 
collapse.  The  failure  event  resulting  from  one  failure  path  can  be  modeled  by  a  parallel 
system. 

Since  there  typically  are  many  different  failure  paths,  each  represented  by  a  parallel 
system,  and  since  failure  can  occur  in  any  one  of  the  failure  paths,  the  entire  system  can 
be  modeled  as  a  giant  series  system  with  parallel  subsystems,  each  representing  a  failure 
path. 

7.3  Structural  Reliability 

In  this  section,  we  consider  terminology  related  structural  reliability,  reliability 
methods,  and  probabilistically  based  structural  design  codes. 
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Reliability 

This  is  the  complement  of  the  peotM^niity  of  fnihse  p^  ie.,  lehi^niity  is  the 
probability  of  survival,  given  by  1  -  p,. 

Safety  Margin 

This  is  the  difference  between  ci^)acity  and  demand,  or  strength  and  load.  Either 
mean  or  characteristic  values  may  be  used  to  determine  the  safety  margin. 


The  basic  concept  of  Level  m  reliability  methods  is  that  a  probability  of  failure  of  a 
structure  always  exists,  and  may  be  calculated  by  integrating  the  joint  probability  density 
function  of  the  variables  involved  ui  the  load  and  strength  aspects  of  the  structure.  The 
domain  oi  integration  is  the  unsafe  regioD  defined  by  the  varbtiiles. 

Because  of  the  difficulties  involved  in  determining  the  joint  density  function  and  in 
calculating  the  multiple  integration.  Level  n  noethods  for  obtaining  the  safety  index  and 
the  related  probability  of  failure  were  introduced.  In  Level  n  methods,  the  probability 
content  of  the  failure  tinmain  is  obtained  using  tqtproximations  to  the  failure  surface. 
F(H(M  and  SORM,  described  elsewhere,  are  Level  n  methods.  Primarily  because  of  the 
iq)proximations  made  to  the  failure  surfaix,  and  also  because  of  iq>proximations  involved 
in  the  inclusion  of  distribution  infotmatitni,  the  probabilities  of  failure  calculated  from 
Level  B  method  ast  not  exact  However,  the  methocte  are  very  effkient  »d  usmdly  a 
good  approximation  is  obtained. 

Level  I  refers  to  safety  factor  based  design  formats  that  are  very  similar  to  traditional 
design  formats  and  safety  check  equations,  except  that  the  safety  factotfs)  are  obtained 
on  the  basis  of  Level  n  methods  to  assure  a  certain  target  reliability  level. 

Safety  Index: 

The  safety  index  is  a  number  that  is  inversely  related  to  the  probability  of  failure. 
The  safety  index  ^  and  the  probability  of  failure  are  related  by 

Pf  *= 

where  is  the  standard  normal  distribution  function.  A  safety  index  of  2.3  translates 
roughly  to  a  probability  of  failure  of  1/100,  3.1  to  1/1000,  and  3.7  to  1/10000.  A  safety 
index  of  zero  corresponds  to  a  probability  of  failure  of  0.5. 


Hasofg-Lind  Safety  Index 

In  the  history  of  structural  reliability  theory,  there  have  been  several  definitions  of  the 
safety  index,  some  fell  from  favor  because  of  a  problem  known  as  lack  of  invariance.  By 
this,  it  is  meant  that  mechanically  different  limit  state  functions  representing  the  same 
physical  fathxro  mode  restated  in  efifferent  values  of  the  safety  index.  The  Hasoler-Lind 
index  does  not  suffer  ftom  the  lack  of  invariance  problem. 

First  Order  Reliability  Method  (FQRMi 

The  essential  steps  in  this  method  of  reliability  aiudysis  for  the  determination  of  the 
probability  of  failure  are: 

a)  The  basic  correlated  random  variables  X  defining  the  limit  state  function  GQO  = 
0,  with  prescribed  probability  distributions,  are  transformed  to  a  set  of 
independent  standard  normal  variaUes  U. 

b)  The  limit  state  surface  gQD  in  the  standard  normal  space  is  approximated  by  its 
tangent  hypeiplane  at  the  point  of  the  limit  surface  closest  to  the  origin.  This 
point  has  the  highest  probability  density,  and  is  called  the  design  point  or  the 
most  probable  failure  point 

c)  The  probability  content  within  the  linearized  failure  domain  is  found  as  an 
estimate  of  the  actual  failure  probability.  The  FORM  probability  of  failure  is 


Pf  = 


where  P  is  the  reliability  index,  which  is  also  the  distance  of  the  design  point  from 
the  origin  in  the  u  space.  The  FORM  reliability  index  is  invariant  for 
mechanically  different  limit  state  functions  representing  the  same  failure  event 

Rackwitz-Fiessler  Transformation 

In  calculating  the  safety  index,  it  is  necessary  to  include  information  related  to  the 
form  of  the  distribution  of  the  basic  variables.  The  tail  of  the  distribution  of  the  random 
variables  is  usually  the  location  where  most  of  the  contribution  to  the  probability  of 
failure  comes  from.  In  the  Rackwitz-Fiessler  transformation,  an  equivalent  normal 
distribution  is  Htted  to  the  tail  of  the  nonnormal  distribution  at  the  most  likely  failure 
point  (design  point).  The  method  requires  the  cumulative  distributions  and  the 
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probability  density  function  of  both  the  actual  distribution  and  the  iKXinal  distributi(»  be 
equal  at  the  design  point 

Second  Order  Rehabilitv  Methods  (SQRhTi 

In  SORM,  the  essential  steps  are  amilar  to  FORM,  except  that  the  limit  state  surface 
in  the  standaxd  normal  ”u"  space  is  approximated  by  a  second  order  tqtproximation  such 
as  a  hypeiparaboloid  fitted  with  its  q>ex  at  the  design  point  Tte  failure  domain 
probability  content  within  the  second  order  iq>proximatic«  is  then  estimated  For 
hyperparaboloids,  the  probability  content  can  be  "exactly”  estimated. 

Safety  Check  Eouation 

In  structural  design,  the  performance  of  the  structure  is  checked  using  safety  check 
equations.  In  the  working  stress  approach  for  fixed  offshore  platforms  as  embodied  in 
API  RP-2A  Recommended  Practice,  fn-  example,  the  maximum  or  yield  strength  is 
divided  by  a  safety  factOT  to  obtain  an  allowable  stress.  Designs  are  then  limited  so  that 
the  tnaaiinum  calculated  stress  under  extreme  operating  loads  does  not  exceed  the 
allowable  value.  This  exanqrle  safety  check  is  of  the  form 

^  ^  D  +  L  +  W  +  otherloadcffects 

where  R  =  nominal  conqionent  strength 
SF  =  safety  factm 

D  =  nominal  gravity  loads  on  conqxments 

L  =  nominal  live  load  effects  on  cmnponents 

W  =  nominal  environmental  load  effects  on  compraents 

Nominal  loads  are  all  combined  with  factors  of  one,  and  constant  safety  factors  1.67  and 
1.25  are  used  for  operating  and  extreme  loadings.  There  are  typically  many  safety  check 
equations  to  be  satisfied  in  a  design,  each  of  which  addresses  a  different  failure  mode  or 
design  concern. 

Partial  Safety  Factor  Format 

A  safety  check  equation  in  a  partial  safety  factor  format  employs  multiple  safety 
factors,  which  may  address  uncertainties  in  component  loads,  resistance,  and  also  failure 
consequences,  non-coincidence  of  peak  loads  from  different  sources,  etc.  Because  there 
is  more  than  one  safety  factor  employed,  the  format  is  more  efficient  in  that  factors  of 
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safety  are  placed  in  a  manner  more  ccnnmensuraie  with  individual  demands  and 
uncertainties.  Also,  the  partial  safety  factors  are  usually  olxained  using  Level  n 
reliability  methods,  consistent  with  a  required  target  reliability  level. 

A  saiiq)le  partial  safety  factor  format  is  that  recommended  in  the  Load  and  Resistance 
Factor  (LRFD)  version  of  API  RP-2A.  This  is  given  by 

>  TfeD  +  T(i,L  +  YwW  +  ... 


where  Rj 


D 

Yd 

L 
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W 

Yw 


s  nominal  strength  or  resistance  of  component  i 
s  partial  resistance  factor  for  component  i 
s  nominal  gravity  or  dead  load  effect 
s  load  factor  for  dead  load 
s  nominal  live  load  effect 
=  load  factor  for  live  load 

=  n(»ninal  environmental  effect  with  prescribed  return  period 
a  load  factor  for  environmental  load 


Each  resistance  factor  4^  is  calculated  as  a  product  of  two  factors,  one  representing 
strength  uncertainty,  and  the  other  taking  into  account  the  consequennce  of  failure  of  the 
component  and  the  structural  system.  The  load  factOTS  7  are  also  calculated  as  a  product 
of  two  factors,  one  representing  uncertainty  in  load  intensity,  and  the  other,  uncertainty 
in  the  related  analysis  ptxx;edures. 

A  partial  safety  factor  format  is  a  Level  I  reliability  based  format  if  the  safety  factors 
employed  are  obtained  from  reliability  analyds  with  a  prescribed  target  reliability. 


Traditionally  in  structural  design,  nominal  or  characteristic  values  are  used  for  the 
basic  design  variables  appearing  in  safety  check  equations.  For  loads,  characteristic 
values  on  the  high  side  of  the  mean  are  typically  used,  while  for  resistance,  characteristic 
values  on  the  low  side  of  the  mean  are  used.  Thus  for  example,  in  ship  design,  safety 
check  equations  involving  yield  strength  use  the  rule  minimum  yield,  which  typically  is 
about  15%  lower  than  the  mean  value.  The  terms  "characteristic"  and  "nominal"  are 
interchangeable,  but  an  occasional  distinction  appears  in  the  literature  where  a 
characteristic  value  refers  to  a  nominal  value  that  is  selected  on  the  basis  of  a  probability. 
For  example,  the  characteristic  yield  strength  may  be  a  S-percentile  value,  i.e.  there  is  a 
95%  chance  that  the  actual  yield  strength  is  greater  than  the  characteristic  value. 
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i^oae  caupranon 

This  is  the  process  of  selecting  a  target  reliability  level  and  a  ctxresponding  set  of 
partial  safety  factors  for  use  in  a  pre^bility  t^ised  design  code.  Reliability  analyses  of 
comparable  past  experience  (existing  structures,  and  systematic  structural  designs  to 
traditional  codes)  are  useful  in  the  code  calilvation  process. 


Code  Optimization 

This  is  the  process  of  selecting  partial  safety  factors  for  use  in  probabilistically  based 
safety  check  equations  in  such  a  manner  that  the  scatter  in  the  reliability  of  structures 
built  to  the  code  is  minimized,  and  centred  around  the  target  value. 


8.  Extrap<4ation  Techniques  for  Design  Loads 


In  this  chapter,  extrapolation  techniqi^s  ftnr  determining  lifetime  extreme  wave  loads 
for  design  axe  identified.  For  purposes  of  discussion,  a.  stnchasric  wave  load,  process  is 
considered.  The  effective  wave  loads  give  rise  to  stress  at  a  point,  which  include  stresses 
arising  from  hull  girder  bending  in  two  planes,  torsion,  external  pressure,  internal  tank 
loads,  etc.  with  proper  accounting  of  phasing. 

Extrapolation  techniques  for  the  wave  load  effect  are  first  considered.  The  definitiem 
of  design  loads  is  subsequently  investigated. 

8.1  Identificatitm  of  Techniques 

There  are  two  broad  classes  of  techniques  for  the  determination  of  the  maximum 
wave  induced  load  over  the  vessel  design  life.  These  are; 

a)  Short  term  techniques,  in  which  the  short  term  statistical  characteristics  of  the 
wave  load  process  in  a  storm  condition  are  used  to  obtain  the  distribution  of  the 
extreme  load,  and  a  characteristic  design  load. 

b)  Long  term  techniques,  in  which  the  long  term  distribution  of  the  wave  induced 
load  is  obtained.  That  distribution  includes  within  it  all  load  peaks  possible 
considering  every  seastate.  A  characteristic  design  load  is  then  defined  based  on 
the  long  term  distribution. 

The  essential  difference  between  the  two  classes  of  methods  is  that  in  the  shon  term 
approach,  the  extreme  load  distribution  in  a  few  high  seastates  is  separately  obtained  for 
each,  and  the  characteristic  design  load  is  typically  taken  as  the  largest  among  values  for 
the  various  seastates,  while  in  the  long  term  approach,  the  design  load  corresponds  to  a 
given  exceedence  probability  (c.g.,  1(H)  on  the  long  term  distribution.  These  two  classes 
of  techniques  are  now  described. 

8.1.1  Short  Term  Wave  Load  Extratxrlation 

If  the  wave  loads  acting  on  a  vessel  can  be  represented  as  a  stationary  Gaussian 
random  process,  which  ir.  usually  an  adequate  a.ssumption  over  the  duration  of  a  sea.state 
lasting  a  few  hours,  then  at  least  two  types  of  methods  are  available  to  predict  the 
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distribution  of  the  maximum  load.  These  two  methods,  among  others,  are  described  in 
detail  by  Mansour  in  [6].  In  the  fust  method,  the  peaks  are  assumed  to  be  statistically 
independent  and  identically  districted,  and  the  distribution  of  the  largest  peak  in  N- 
peaks  is  determined  using  classical  order  statistics.  In  the  second,  conventional 
upcrossing  analysis  is  used  for  determining  the  extreme  value  distribution. 

A.  Distribution  of  largest  peak  by  OTder  statistics 

The  distribution  of  the  largest  peak  in  a  sequence  of  N-peaks  can  be  determined  using 
standard  order  statistics.  Consider  a  sequence  of  random  variables,  Zj,  Zj,  ...  z^ 
representing  the  peaks  of  a  load  on  a  marine  structure.  Assuming  that  these  peaks  are 
identically  distributed  and  statistically  independent,  the  cumulative  distribution  function 
of  the  largest  one  is  given  by 

(z)  =  P[max  (z,,  Zj, ...  z^)  ^  z] 

=  [F.  (z,e)]N 

where  F,  (z,e)  is  the  initial  comulative  distribution  function  of  the  load  peak  (maxima) 
and  e  is  the  spectral  bandwidth  parameter  defined  from 


mo®4 

m„  =  05“  S(05)do5;  n  =  0, 2, 4 

Here,  0)  is  the  radian  frequency.  The  probability  density  function  (pdf)  of  the  largest 
peak  is  determined  by  differentiating  the  c.d.f.  with  respect  to  z,  thus 

f^  (z)  =  N[F,  (z,e)]N->  .  f,  (z,E) 

where  4  (z,e)  is  the  initial  p.d.f.  of  the  load  peaks.  For  an  arbitrary  bandwidth  process, 
the  initial  distribution  of  peaks  within  a  short  term  seastate,  considering  positive  maxima 
alone,  has  been  derived  by  Ochi  (J.  Ship  Research,  1973).  For  a  defrnition  of  positive 
and  negative  maxima  and  positive  and  negative  minima,  see  Figure  8.1.  In  the  narrow 
banded  case,  the  conventional  Rayleigh  density  and  distribution  functions  apply. 
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Based  on  the  Ochi  distribudon  and  order  statistics,  it  can  be  shown  that  the  modal 
value,  Le.,  the  most  probable  maximum  load  in  N-peaks  is  aiyroximated  by 


The  i^tproximatiog  was  derived  by  Octe  coasiderifig  hffge  N  and  6  £  (X9.  It  can  be 
shown  that  there  is  a  63%  chance  that  the  largest  lespmise  will  exceed  the  modal  value. 
Other  percentile  values  of  the  extreme  value  distributimi  were  also  obtained  by  Ochi,  in 
terms  of  a  "risk  parameter"  a.  He  chooses  a  very  small  number,  a  (e.g.,  0.01)  and 

A 

obtains  a  non-dimensional  extreme  value  4n  ^ 


Extreme  value  of  maxima  £ 

inNpeaks 

For  e  £  0.9,  N  lai:ge,  and  a  small,  it  can  be  shown  that 

A 


The  dimensional  extreme  value  is  equal  to  the  non-dimensional  extreme  value  multiplied 

byVi% 


B.  Extreme  value  distribution  based  on  upcrossings 

The  distribution  of  the  largest  peak  can  be  determined  from  upcrossing  analysis  of  a 
rime  history  of  a  stationary  random  process  instead  of  the  peak  analysis  described  above. 
Principles  behind  the  upcrossing  analysis  are  described  by  Mansour  (Ship  Structure 
Committee  Report  351)  and  will  not  be  repeated  here.  The  essential  problem  is  one  of 
determining  the  first  passage  of  a  random  process  x(t)  of  a  level  "a"  within  a  given  time 
interval  T.  Based  on  a  level  oossing  analysis,  assunung  that  the  individual  level  arrivals 
are  independent  and  Poisson  distributed,  it  can  be  shown  that  the  cumulative  distribution 
function  of  the  largest  x  value,  denoted  Z,  is 

Fz  (a)  =  exp  (-v,  (a)  T) 
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where  v,  (a)  is  the  expected  number  of  level  crossinp  per  unit  time.  This  is  givoi  by 

v.d)  .  v,exp(-5^j 

In  the  above,  Vg  is  the  zero  crossing  rate,  which  for  a  nanow  band  process,  is 

V 

®  2ji  \mg 

The  above  cumulative  distiibutitMi  function  for  die  largest  value  ignores  the  tendency  fOT 
upcrossings  to  occur  in  clumps,  because  of  the  assumption  of  indqiendence.  The 
solution  overpredicts  extreme  values.  To  consider  clurxqnng,  an  upper  bouiul  oivelc^  to 
the  given  process  can  be  constructed,  and  the  first  passage  probability  for  the  envelc^ 
process  obtained.  The  upoossing  rate  V|t(a)  for  the  envelope  of  a  Gaussian  process  is 
given  in  standard  structural  reliability  textbodks  as 

Vj(a)»V2x  I — - Vx(a) 

\  rngm,  .y/mg 

In  general,  this  upcrossing  rate  will  not  lead  to  a  decreased  bound,  since  the  envelope 
may  have  excursions  above  the  level  without  there  being  actual  process  upcrossings. 
Such  crossings  are  termed  "empty",  while  otherwise  they  are  called  "qualified" 
upcrossings,  a  terminology  devised  by  Vanmarcke  (ASME,  J.  Applied  Mechanics,  March 
1975).  Vanmarcke  obtained  an  estitruue  of  qualified  excursions,  which  was  later  refined 
using  a  Slepian  regression  model  by  Didevsen  and  Lindgren  (J.  Sound  and  Vibration, 
1988). 

To  date,  the  Didevsen  and  Lindgren  solution  is  the  best  available.  Based  on  it,  the 
cumulative  distribution  function  of  the  maximum  value  for  an  ergodic  Gaussian  narrow 
band  wave  load  process,  becomes.  (Cramer  and  Friis  Hansen,  "Stochastic  Modeling  of  the. 
Long  Term  Wave  Induced  Responses  of  Ship  Structures,"  submitted  to  Journal  of  Marine 
Structures): 
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exp\ 


Tr,(a)v,(a) 


wtoe  (a)  was  inevioosly  defined,  V  is  the  level  value,  and  ty  (a)  is  given,  for 
moderate  spectral  skewness,  fiom 


where 


l-ry(a)  =  2r0(ti) 


<b  Yijc - - - 
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The  extreme  value  analysb  based  oa  the  upcrossing  rate,  as  d>tained  above,  provides 
a  cumulative  distribution  function  of  the  extreme  value,  accounting  for  clumping  of 
peaks.  It  is  derived  for  a  narrow  band  ergodic  Gaussian  wave  load  process,  although 
based  on  simulation  coirparisons,  it  seems  ^licable  to  relatively  wide  band  processes 
also.  It  is  worth  stating  that  the  probability  doisitv  function  of  the  maximum  value  has 
not  been  obtained. 


C.  Calculation  of  the  short  term  extreme  values 

Short  term  extreme  values  based  on  the  peak  or  level  crossing  analyses  are  calculated 
seastate  by  seastate  for  several  extreme  wave  conditions.  Within  a  seastate,  the  extreme 
values  depend  on  (are  conditional  on)  vessel  heading  and  speed.  Typically  in  treating 
low  frequency  wave  induced  loads,  the  speed  within  a  seastate  is  assumed  constant  and 
extreme  values  conditional  on  difierent  wave  headings  are  obtained.  The  extreme  value 
for  the  seastate  is  obtained  by  unconditioning  with  respect  to  vessel  headings,  i.e.,  the 
wave  load  extreme  values  for  each  heading  are  multiplied  by  the  heading  probabilities 
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and  added  The  largest  characteristic  extreme  load  among  all  seastates  cauidered  may 
be  used  as  a  design  load 


In  the  long  term  tqtpioach  to  the  entire  density  or  distribution  lunctkm  the  wave 
induced  load  is  obtained,  oxtsidenng  the  following: 

(i)  Frequency  of  occurrence  of  various  sea  states. 

(ii)  Frequency  of  occurrence  of  various  spectral  shapes  witiun  each  sea  addition. 

(iii)  Ship  route  and  frequency  of  encountering  each  seastate  and  qrectral  shape. 

(iv)  Frequency  of  occurrence  of  various  vessel  headings. 

(v)  Frequency  of  occurrence  of  various  vessel  speeds. 

(vi)  Frequency  of  occurrence  of  various  ship  loading  conditicms. 

(vii)  The  expected  number  of  load  cycles  for  a  given  sea,  wave  spectral  shape, 
speed  and  heading. 

The  consideration  of  various  spectral  shapes  within  a  seastate  is  characteristic  of  some 
procedures  based  on  seastate  groups,  where  the  seastates  possible  in  the  long  term  are 
grouped  into  a  small  number  of  "weather  groups”.  An  exampb  will  be  given  later. 

Taking  the  various  factors  noted  into  consideration,  the  probability  density  function 
of  the  load  peaks  applicable  to  the  long  term  response  can  be  written  for  each  ship 
loading  condition  as: 

ZZZZ  “.PiPjPkP/f.W 

f  (X)  =  - 

"‘PiPiPiP/ 

i  j  k  < 

where  f^(x)  is  the  probability  density  function  for  the  load  peaks  in  the  short  term,  and 
n^  is  the  associated  number  of  peaxs  per  unit  time.  For  a  narrow  band  process,  n^  is 
obtained  based  on  the  Rayleigh  density  for  peaks  in  the  short  term,  as 
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n.=f 

The  weighting  factor  pj  represents  the  expected  occurrence  probability  for  the  sea 
;x)odition,  pj  for  the  wave  specinnn  shqre,  i\  ftH*  headings  in  waves  in  a  given  sea  and 
spectrum  shtqre,  and  p^  for  speed  in  a  given  sea«  spectrum  shape  and  heading.  The  total 
number  of  responses  expected  during  the  vessel  life  then  becomes 


NT  =  XXZZ^"*PiPiPk  P/)xTx60* 

i  j  k  < 

where  T  is  the  total  sea  exposure  time  in  hours.  The  formula  for  the  probability  density 
function  and  the  total  number  of  cycles  a  plies  to  wide  band  short  term  processes  also, 
with  n^  and  f^(x)  appropriately  calculated.  The  cumulative  distributitm  function  of  the 
wave  load  in  the  long  term  is  also  similarly  obtained. 

It  is  worth  reiterating  that  in  the  long  term  approach,  distribution  and  density 
functions  in  the  long  term  are  obtained  by  weighting  and  adding  the  short  term  density 
and  distribution  functions.  The  short  term  density  and  distribution  functions 
corresponding  to  the  peaks  (e.g.,  Rayleigh  distributitm)  are  generally  used.  For  the  long 
term  distribution  thus  obtained,  the  probability  scale  includes  each  peak  or  load  cycle. 
The  load  corresponding  to  a  l/N-i-  exceedence  level  is  often  used  as  the  design  load.  If 
Nf  =  10*,  as  is  the  case  in  merchant  ships,  the  exceedence  level  is  lO-*,  and  the 
corresponding  "10-*  load”  is  used  as  a  design  load. 

The  Weather  Group  Approach 

In  the  typical  long  term  approach,  a  wave  scatter  diagram  for  the  long  term  is  used. 
Each  bin  in  the  scatter  diagram  characterizes  a  seastate  defined  by  a  significant  wave 
height,  a  spectral  period,  and  an  associated  occurrence  probability.  In  calculating  the 
wave  loads,  one  analytical  seaspectrum  such  as  that  due  to  Bretschneider  or  ISSC,  is  used 
for  each  bin  of  the  scatter  diagram. 

In  an  alternate  approach,  the  long  term  wave  environment  is  discretized  into  weather 
groups,  with  associated  probabilities.  For  the  average  North  Atlantic,  Lewis  in  1967 
suggested  the  following  weather  groups  and  associated  fiequencies  of  occurrence; 
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H,n.  feet 

%  occurrence 

10 

84.54 

20 

13.30 

30 

2.01 

40 

0.14 

48.2 

0.01 

In  each  weather  group,  more  than  one  preselected  wave  spectrum  (typically  about  10) 
must  be  used  for  the  short  teim  wave  load  calculations.  The  spectral  forms  used  are 
typically  based  on  measurements,  and  represent  a  range  of  wave  peak  frequencies.  The 
long  term  distribution  is  constructed  from  the  short  term  distributions.  In  the  process, 
some  weather  group  methods  may  assunw  each  spectral  form  within  a  weather  group  to 
have  predefined  probabilities  of  occurrence.  Others  may  use  additional  (predefined) 
information  on  the  spread  of  short  term  mean  square  values  within  a  weather  group. 

The  weather  group  approach  may  also  be  termed  ..  "spectral  family"  !q>proach. 
Spectral  families  for  the  North  Atlantic,  which  is  the  design  wave  environment  foi 
merchantships,  have  also  been  provided  by  Ochi,  SNAME  Transactions,  1978.  A 
weatheigroup  approach  based  on  wave  ^)ectral  measurements  in  the  North  Atlantic  is 
used  by  the  American  Bureau  of  Shipping  for  vessel  smictural  assessment  for 
unrestricted  service. 
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Fig.  8.1  Explanatory  Sketch  of  a  Random  Process 
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8.2  Detemiinatioii  of  Design  Loads 


Methods  for  the  extrapolation  of  wave  induced  load  were  considered  in  the  {Rtvious 
section.  In  this  section  we  consider  how  design  loads  are  defined.  There  e«:esdslly  are 
two  possible  criteria  for  the  definition  of  design  loads.  These  are  to 

a)  Select  the  loads  such  that  a  certain  level  of  exceedence  is  acceptable  on  the  basis 
of  either  short  or  long  term  procedures. 

b)  Select  the  loads  such  that  the  structural  reliability  level  considering  one  or  more 
limit  states  is  acceptable. 

We  illustrate  the  two  procedures  considering  a  Stillwater  load,  a  wave  induced  load  and  a 
strength  variable.  The  problem  of  treating  combined  loads  for  the  same  purpose  of 
identifying  design  loads  is  an  advanced  one,  and  is  in  fact  part  of  a  ship  structure 
committee  research  project  cm  Load  Combinations,  SR-1337.  Our  more  basic  treatment 
considers  the  Stillwater  load,  wave  load  and  strength  to  be  independent  of  one  another. 

8.2.1  Selection  of  Maximum  Load  Effect  for  Design 

With  a  single  wave  load  present,  there  is  a  one  to  one  correspondence  between  the 
load  and  the  load  effect  In  this  context  the  Stillwater  load  is  not  specifically  considered. 
Because  it  is  essentially  constant  over  voyages  that  last  days  or  a  noonth,  its  inclusion  or 
consideration  does  not  pose  a  difficult  problem.  The  only  question  to  be  answered,  then, 
is  how  to  determine  the  maximum  expected  wave  load  in  the  lifetime  of  the  vessel.  Such 
load  is  pertinent  to  structural  design  for  extreme  loads. 

We  previously  described  two  methods  for  obtaining  the  distribution  of  the  largest 
wave  load  peak,  either  by  using  order  statistics  or  by  level  crossing  analysis.  These  two 
classes  of  methods  apply  to  a  short  term,  i.e.,  seastate  by  seastate  analysis.  We  also 
described  methods  for  the  construction  of  the  long  term  wave  load  distribution, 
considering  every  load  peak  in  each  seastate.  The  following  are  the  typical  ways  of 
defining  the  extreme  wave  load  for  design,  based  on  the  above  approaches: 

Short  Term  Analysis: 

In  design,  the  largest  wave  load  is  defined  considering  the  most  probable  value  of  the 
wave  load  distributions  in  each  possible  seastate.  The  number  of  short  tnm  wave  load 
peaks  N  is  computed  from  the  zero  crossing  period  for  each  seastate.  The  design  wave 
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load  is  the  largest  among  the  set  of  short  term  most  probable  extrenw  wave  loads  far  the 
selected  seastates. 

The  seastates  should  be  selected  on  the  basis  of  an  acceptable  return  period  and/or 
acceptable  probability  of  the  ship  encountering  such  seastates.  The  latter  depends  cm  the 
operational  life  and  the  route  of  the  ship.  Reference  [6]  describes  tecfaniqnes  for 
computing  probability  of  encountering  a  seastate  of  a  specified  return  period,  as  well  as 
techniques  for  determining  a  seastate  with  a  specified  return  period  based  on  wave  data. 

Long  Term  Analysis 

In  this  method,  the  design  value  is  taken  to  be  the  largest  wave  load  with  an 
exceedence  probability  of  1/N,  N  being  the  total  number  of  wave  load  peaks.  In 
calculating  N,  and  in  obtaining  the  long  term  distribution,  each  wave  load  peak  possible 
is  considered.  If  the  total  number  of  load  peaks  is  iO^  in  20  years,  for  example,  the 
design  value  is  the  10**  exceedence  level  value  from  the  long  term  distribution.  This 
value  is  said  to  occur  once  in  the  lifetime  of  the  vessel. 

While  not  usual,  risk  parameter  can  also  be  included  in  the  long  term  approach.  The 
design  value  of  the  wave  load,  ,  is  then  determined  such  that 

l-F(ZN)=f 

where  a  is  the  risk  parameter,  e.g.,  0.01,  N  is  the  total  number  of  cycles  (i.e.,  wave  load 
peaks)  in  the  long  term,  and  F(Zf()  is  the  cumulative  distribution  function  of  the  long 
term  wave  load. 

8.2.2  Design  for  a  Target  Reliabilitv  Level 

Probabilistic  methods  provide  a  mechanism  for  obtaining  extreme  design  loads  for  a 
structure  with  the  required  target  reliability  or  failure  probability.  The  design  safety 
check  equation  for  the  limit  state  may  take  the  conceptual  form 

<>  C^y.D.  +  y,D, 

where  ()>  is  the  strength  partial  safety  factor,  and  7,  and  7^  arc  the  still  water  and  wave 
load  partial  safety  factors.  The  C,  D,  and  are  characteristic  values  of  the  strength, 
still  water  and  wave  loads.  The  seastate  that  defines  D,  was  previously  identified.  The 
problem  is  then  one  of  determining  ({>,  7,  and  considering  the  uncertainties  in  strength 
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and  loads,  such  that  a  target  reliability  level  is  achieved.  Level  1  reliability  methods  can 
be  used  in  this  process.  The  derivation  of  the  partial  safety  factors  associated  widi  each 
design  variable,  including  the  loads,  fra*  a  target  reliability  level  is  described  in  Part  1  of 
this  report.  For  additional  discussion  of  such  procedures,  the  reader  is  referred  to 
Mansour  [6]. 
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9.  Serviceability  Limitstatcs 


This  chapter  pertains  to  ufendfication  and  descr^tion  of  hnpottant  serviceabifity  Emit 
states.  By  definition,  a  serviceabUity  Emit  state  is  associated  with  constraints  on  die 
structure  in  terms  of  requirements  such  as  maximum  deflection  of  a  member,  critical 
buckling  loads  that  cause  elastic  buckling  of  a  plate  element,  or  local  cracking  due  to 
fatigue.  The  limit  state  manifestatimis  are  typicaEy  of  aesthetic,  fiincticHial  or 
maintenance  concern,  but  do  not  normally  lead  to  overall  collapse.  The  following 
serviceabUity  limit  states  are  now  considered. 

(a)  serviceabUity  limit  state  associated  with  critical  buckling  stresses 

(b)  serviceability  limit  stale  associated  with,  fatigue 

9.1  ServiceabUity  Limit  State  for  Plate  Budding 

Plate  elements  in  a  ship  huU,  such  as  between  Itmgitudinals,  can  buckle  under  appUed 
loads  in  either  the  linear  elastic  or  inelastic  range  of  material  behavior.  A  plate  that 
buckles  in  the  linear  elastic  regime  wUl  essentially  regain  its  original  configuration  when 
unloaded.  On  the  other  hand,  a  plate  that  buckles  in  the  inelastic  regime  may  suffer  some 
permanent  set  upon  unloading.  The  appUed  stress  that  defines  the  lower  limit  of  the 
inelastic  regime  is  that  corresponding  to  the  material  proportional  Emit  Thus  the  so- 
caUed  inelastic  regime  includes  nonlinear  elastic  and  plastic  behavior. 

Buckling  of  plate  elements  in  the  linear  elastic  regime  is  generaUy  acceptable  in 
longitudinaUy  framed  vessel  huUs,  although  it  is  rare  that  the  designer  intentionaUy 
designs  the  structure  to  behave  so.  The  major  exceptions  to  this  occur  in  passenger 
vessels  and  car  carriery  where  the  plating  on  decks  above  tire  weather  deck  from  stress 
considerations  alone  can  be  relatively  thin,  their  main  function  being  to  provide  the 
required  weather  and  water-tightness.  In  such  cases,  it  is  efficient  for  the  designer  to 
allow  linear  elastic  plate  buckling  to  occur,  the  result  being  a  tighter  structure  than  would 
otherwise  be  the  case,  and  also  less  topside  weight. 

Depending  on  the  philosophy  of  the  profession  and  the  organizations,  buckling  of 
plate  elements  in  the  inelastic  regime  may  or  may  not  be  allowed,  the  primary 
consideration  being  aesthetic.  From  a  material  utilization  point  of  view,  plate  thicknesses 
can  be  reduced  if  an  amount  of  permanent  set  is  allowed. 

In  discussing  serviceability  limit  states  involving  plate  buckling  in  longitudinally 
framed  vessels,  the  following  nomenclature  is  adopted:  The  plate  long  dimension 
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Qength)  is  assunKd  to  be  parallel  to  the  x  axis,  or  the  vessel  longitudinal  diiecticm,  and  is 
labeled  "a".  The  plate  width  ot  small  dimensitHi  is  taken  parallel  to  the  y  axis  <x  vessel 
transverse  direction,  and  is  labeled  "b".  The  plate  aspect  ratio  aA>  is  always  larger  than  or 
equal  to  unity.  The  plate  thickness  is  denoted  ”t".  The  plate  element  is  ccmsidered  under 
uniform  inplane  conqnession,  either  in  the  longitudinal  directicm  (the  so-called  long  plate 
case)  or  in  the  transverse  (hiection  (the  so-called  wide  {due  case).  Another  load  case 
ccmsidered  is  the  plate  under  uniform  edge  shear.  The  serviceability  limit  state  is  reached 
when  the  applied  stress  equals  o,  or  Op,  where  the  limit  o,  applies  in  the  linear  elastic 
range,  and  Op  applies  in  the  inelastic  range. 

Uniform  Compression 

For  long  plate  compression, 

Jc"E 

“■  12(I-v=) 

where  k  =  4  for  simply  supported  edges.  For  other  edge  conditions,  the  buckling 
coefficient  K  can  be  obtained  from  the  attached  Figure  9.1.  If  ^  Gp^^,  the 
proportional  limit. 


Otherwise, 


0lL(<yY-<yil)  +  0(*' 


In  the  above,  Gy  is  the  material  yield  strength. 
For  wide  plate  compression. 


I,  "'E  fiY 


where  k 


for  simply  supported  edges.  If  Gcr  S  GpL,  the  following  applies. 


• 
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Odwrwise, 


o-*o. 


Oil  (gy -Oil) 


Edge  Shear 

The  critical  buckling  stress  is  given 

-V 

*^12(l-v*)lbJ 

where  k  =  5.34+4^— j 

for  sinq>ly  supported  edges.  If  the  edges  can  be  considered  clanged,  dte  buckling 
coefficient  k  takes  the  form 

k*8.9«+5*6^^J 


In  the  Unear  elastic  range,  that  is,  if  ^  <W^> 

ge  “  ‘*CR 


Otherwise,  the  limit  stress  is 


_ 3T,tat^ 

^'"gH.(gy~gH.)+3ta,* 

where  ty  is  the  shear  yield  stress,  equal  to  Gyhfs. 

The  above  solutions  defining  the  serviceability  limit  states  under  uniform  edge 
compression  or  shear  are  based  on  classical  buckling  theory.  Further  reference  is  made 
to  [8].  The  limit  stresses  beyond  the  proportional  limit  are  based  on  tangent  modulus 
corrections  due  to  Bleich.  The  interested  reader  is  referred  to  Bleich's  book  on  "Buckling 
Strength  of  Metal  Structures",  published  by  McGraw  Hill,  1952.  With  a  tangent  modulus 
correction  11  included,  the  limit  stress  can  be  written  in  the  following  form; 


where  11 »  f(E|3).  Et  being  the  tangent  modulus.  The  functional  relatkmships  defixung  ii 
are  different  for  the  long  plate,  wide  plate  and  shear  cases.  Hrace  the  coneqmnding 
are  also  of  different  fmns. 

9.2  Serviceability  Limit  State  for  Fat^e 

The  fatigue  limit  state  is  associated  with  the  danuiging  effect  of  repeated  loading 
which  may  lead  to  loss  of  a  specific  fiiiKnion,  maintenance  costs,  and  in  certain  cases  to 
ultimate  collapse.  That  fatigue  cracks  in  ships  are  more  a  maintenance  than  a  safety 
concern  is  essentiaSy  due  to  the  ductifity  of  ship  steels.  Fatigue  cracks  do  occur  in 
complex  structures,  and  design  against  fatigue  (i.e.,  procedures  to  limit  fatigue  cracking 
to  acceptable  levels)  is  important 

There  are  various  possible  ways  of  computing  the  fatigue  damage  in  a  vessel  subject 
to  a  specified  long  term  wave  environment  According  to  [6],  the  different  methods  may 
be  classified  as  those  based  on 

(a)  wave  height  history 

(b)  stress  nuige  history 

(c)  the  entire  scatter  diagram 

This  method  of  classification,  further  explained  in  the  attached  Figure  9.2,  is  based  on  the 
level  of  detail  in  the  treatment  of  the  environment  Other  types  of  classification  axe  also 
possible,  e.g.,  S-N  curve  based  methods  as  q>posed  to  fracture  mechanics  based  methods, 
design  stage  methods  in  contrast  to  design  chedting  mediods,  and  so  on.  The 
formulation  of  the  fatigue  limit  state  will  depend  on  the  details  of  the  method  used.  In 
this  section,  the  formulation  described  is  the  one  used  in  section  4.2.3  of  this  report 

The  limit  state  formulation  is  based  on  S-N  curves,  which  describe  the  number  of 
constant  amplitude  stress  cycles  to  failure,  as  a  function  of  the  fluctuating  stress 
amplitude.  The  curve  is  written  in  the  form 

NAS“  =  C 

where  N  is  the  number  of  cycles  to  failure,  AS  is  the  constant  amplitude  stress  range,  and 
m  and  C  are  slope  and  intercept  related  constants.  For  design  purposes,  C  is  chosen  so 
that  the  S-N  cmve  forms  a  "lower  bound"  to  die  experimental  data.  One  typical 
statistical  way  of  defining  C  is 
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where  C  conespoods  to  the  median  S-N  curve,  and  is  the  standard  deviatioD  of 
N.  Each  generic  structural  detail  type  has  an  S-N  curve.  For  a  collection  of  S-N  curves 
typical  of  ship  structural  detail  situations,  the  reader  is  referred  to  Munse's  Ship  Structure 
Committee  report  SSC-318,  Tadgue  Characterizati(»  of  Fabricated  Details  for  Design". 

The  wave  environment  is  described  oxx^letely  by  the  set  of  seastates  and  their 
pfobatnHties  of  occunence  as  defmed  m  a  scatter  (ha^un.  For  each  seastaie,  the  seess 
distribution  can  be  considered  Rayleigh  distributed,  assuming  that  the  wave  induct 
stress  process  is  narrow  band  and  zero  mean  Gaussian.  The  Raylei^  density  is  of  the 
form 

where  is  the  zero  moment  of  the  stress  spectrum  in  seastaie  "j”.  This  nooment  is  also 
equal  to  the  mean  square  value  of  the  stress  process.  The  zero-crossing  frequmcy  of  the 
stress  process  in  hertz  (cycles/secmid)  is  given  by 


where  Xjj  is  the  second  moment  of  the  stress  spectrum  for  the  seastate.  If  the  time  spent 
in  the  seastate  is  Tpj,  where  T  is  the  total  time  period  and  pj  is  the  probability  of 
occurrence  of  the  seastate,  the  number  of  stress  cycles  associated  with  the  seastate  is 

Also,  the  number  of  cycles  associated  with  a  stress  interval  ds  is  [f,(s)ds]  •  Tpj  fj. 

The  fatigue  damage  associated  with  the  seastate  "j"  can  then  be  calculated  using  the 
Miner  linear  cumulative  damage  hypothesis.  The  damage  is  given  by 


where  N|(AS)  is  the  number  of  cycles  to  failure  at  the  q)ecified  stress  range  AS  as 
detem^ined  firam  the  S-N  oisve.  Suhsdturing  for  N((AS)ir  the  abow  may  be 

rewriten  as  follows: 


D 

Here,  the  integral  has  been  evaluated  by  substituting  the  Rayleigh  density  for  f(s).  From 
this,  and  upon  substituting  for  fj,  the  total  damage  in  time  T,  for  j  seastates,  may  be 
obtained  as 


-  ^  Pj  2“  f 
j  r  J 


s“  f(s)ds 


The  above  equation  defines  the  fadgue  damage  from  the  entire  scatter  diagram,  for 
the  time  period  T.  If  the  Palmgren-Miner  damage  sum  at  failure  is  denoted  Af,  the  time 
to  failure  may  then  be  obtained: 


Equation  4.5  of  the  text  is  directly  obtainable  from  the  above  equation  for  time  to  failure. 
That  equation  also  includes  a  stress  inaccuracy  term  B  which  represents  the  "modeling 
error"  in  the  procedures  used  to  compute  the  wave  induced  stress. 

The  above  definition  of  the  fatigue  limit  state  equation  in  terms  of  the  time  to  failure 
assumes  that  the  stress  process  within  any  seastate  in  the  scatter  diagram  is  narrow 
banded.  A  correction  for  the  possible  wide  banded  nature  of  the  process  is  available,  see 
Wirsching  and  Light,  ASCE  Journal  of  Structural  Division,  Vol.  106,  No.  ST7,  July 
1980.  The  wide  band  correction  was  derived  by  Wirsching  and  Light  using  rainflow 
counting  on  simulated  time  histories  of  differing  bandwidths  to  obtain  the  stress  range 
histogram  and  then  computing  the  fatigue  damage,  which  then  was  compared  to  that 
predicted  from  the  narrow  band  assumption.  The  importance  of  the  refinement  obtained 
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by  innlu<<ing  the  correcticm  is  relatively  small  whm  ccxnpared  to  the  inaccuracies 
introduced  by  the  stress  modeling  error  in  particular.  Also,  the  cctrection  assumes  the 
«»5tim«te!t  obtained  by  a  rainflow  count  based  pjrocedure  to  be  die  correct  ones. 
Nevertheless,  the  rainflow  correction  pnovides  a  means  for  obtaining  a  fatigue  damage 
that  is  scNxiewhat  more  realistic  than  that  calculated  using  the  narrow  band 
assumptions.  For  typical  stress  time  histories  in  ship>s,  the  effect  of  the  oocrecticNi  is  to 
reduce  the  calculated  damage. 
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i  :  3  4  s 
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Fig.  9.1  Buckling  Q)efficients  for  Plates  in  Uniaxial  Compression 
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Fig.  9.2  Qassificadon  of  Fbur  Basic  Methods  of  Computing  Fatigue  Damage 


10.  Limit  States  Associated  with  Lifetime  Extreme  Loads 


The  aim  of  this  chapter  is  to  identify  and  describe  the  appropriate  limit  states 
associated  with  lifetime  design  extreme  loads.  The  following  global  limit  states  are 
considered: 

(a)  Hull  girder  initial  yield  limit  state 

(b)  Hull  girder  fully  plastic  limit  state 

(c)  Hull  girder  collapse  limit  state 

The  strength  associated  with  the  following  local  buckling  limit  states  are  also  considered: 

•  Column  and  beam  column  buckling  of  longitudinals 

•  Torsional/flexural  buckling  (tripping)  of  longimdinals 

•  Grillage  buckling  of  longitudinals  together  with  transverse  beams 

The  global  limit  states  apply  to  the  hull  girder  as  a  whole.  The  local  limit  states 
apply  to  portions  of  the  hull  girder,  e.g.,  longitudinals  between  transverses,  longitudinals 
and  associated  flange  plating  between  transverses,  or  gross  panels  consisting  of 
longitudinals  and  transverses.  Plate  bttckling  per  se  is  not  considered,  except  to  the 
extent  it  reduces  the  effective  flange  plate  acting  together  with  the  longitudinals. 

Global  and  local  behavior  are  interlinked,  and  an  argument  may  be  made  that 
consideration  of  global  behavior  alone  is  sufficient  provided  the  consideration  is  detailed 
enough.  Nevertheless,  a  two  level  approach  is  used  because 

(a)  Separate  consideration  of  local  behavior  affords  the  designer  more  control  over 
material  deployment. 

(b)  Local  behavior  is  often  indicative  of  global  behavior. 

(c)  A  two  level  limit  state  design  procedure  is  more  consistent  with  present 
conventional  design  practice. 

(d)  The  fact  that  local  behavior  has  been  controlled  in  design  to  acceptable  levels  can 
lead  to  procedural  simplifications  in  the  consideration  of  global  behavior.  A 
simple  example  is  a  situation  where,  if  buckling  cannot  occur  in  any  local  portion 
of  a  longitudinally  framed  tanker  to  a  given  load  level,  global  considerations  can 
usually  exclude  buckling,  again  up  to  that  load  level. 
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10.1  Hull  Girder  Limit  States 


10.1.1  Initial  Yield  Limit  State 

In  this  limit  state,  hull  girder  behavior  as  a  beam  is  considered.  The  geometric 
property  that  characterizes  hull  girder  behavior  is  its  section  modulus.  It  is  assumed  that 
under  the  applied  extreme  bending  moment,  the  various  elements  of  the  hull  cross  section 
remain  stable,  i.e.,  no  buckling  occurs.  Tlie  stress  at  any  location  'y'  above  the  neutral 
axis  of  the  hull  girder  (see  Figure  10.1)  cross  section  is  given  by 

„  M(x)y 
*  I(x) 

uteie  ;  tfaepnmaiyloiigitiidinalbeitdi^sticsattkKaikmx 

y  :  distance  from  neutral  axis  of  section  to  the  location  where  the  stress  is 
computed 

M(x) :  External  bending  moment  at  longitudinal  location  x 

I(x)  :  moment  of  inertia  of  the  cross  section  at  longituouud  location  x 

Note  that  I(x)/y  is  the  elastic  section  modulus,  and  the  stress  is  maximum  for  minimum 
I/y,  i.e.,  maximum  'y'  distance.  One  can  define  the  first  yield  moment  for  the  cross 
sectkm  as  fellows  (location  parameter  'x*  ormtted): 

=SM,a, 

where  Mfy  is  the  first  yield  moment,  SM^  is  the  minimum  elastic  section  modulus  at  the 
location  of  maximum  bending  moment,  and  is  the  material  yield  strength.  This 
expression  assumes  elastic  behavior  until  the  stress  at  the  extreme  Hbers  reach  yield.  The 
first  yield  moment  is  in  principle  different  for  different  longitudinal  locations.  At  any 
location,  the  first  yield  moment  is  only  realized  if  buckling  does  not  occur.  Nevertheless, 
the  first  yield  limit  state  is  commonly  used  as  a  convenient  strength  characterization 
parameter  in  ship  hull  design. 

10.1.2  Fully  Plastic  Lirmt  State 

In  the  first  yield  limit  state,  the  limit  strength  was  defined  as 
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Mfy  =  SM^a, 


where  SM,  is  the  minimum  elastic  section  modulus,  usually  given  at  any  cross  section  as 
where  I  is  the  moment  of  inertia  of  the  cross  section  and  'y'  is  the  distance  from  the 
neutral  axis  to  the  extreme  fiber  (deck  or  bottom).  The  stress  distribution  is  linear  from 
the  neutral  axis  to  the  location  under  consideration,  and  only  the  Tnaximiiin  stress  at  the 
extreme  fiber  is  at  yield. 

In  contrast,  in  the  case  of  the  fully  plastic  limit  state,  the  entire  cross  section  of  the 
hull  including  sides  has  reached  yield.  The  changes  in  stress  distribution  from  the  first 
yield  to  die  fully  plastic  limit  state  are  sketched  in  Figure  1(12  for  an  idealized  box  girder 
cross  section.  The  following  are  assumed; 

a)  Elastic  perfectly  plastic  material  behavior 

b)  No  buckling 

c)  The  applied  external  moment  does  not  change  direction 

For  the  box  girder  cross  section,  the  fully  plastic  moment,  defined  as  the  internal 
resisting  nooment  with  the  entire  cross  section  at  yield,  may  be  written  as 

Mp  S'  CySMp 

where  oy  is  the  material  yield  strength,  and  SMp  is  a  plastic  section  modulus.  It  can  be 
shown  that 


SM,  =  A„  g+ A,(D-g)+2A,^|-g+^j 

where  =  cross  sectional  area  of  one  hull  side,  the  thickness  being  t,^ 

Ab  =  cross  sectional  area  of  bottom 
A])  =  cross  sectional  area  of  deck 
D  =  depth 

The  areas  include  stiffening  and  plating.  The  variable  'g'  represents  the  distance  from  the 
center  of  deck  area  to  the  plastic  neutral  axis.  The  plastic  neutral  axis  is  defined  by  a 
condition  that  the  areas  above  and  below  it  are  equal,  for  purposes  of  force  equilibrium. 
The  location  of  plastic  neutral  axis  is  defined  by 
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D  AA, 

For  more  compticaied  cross  sections  anc^  if  more  than  one  material  is  osed  in  die 
hull  cross  section,  the  fully  plastic  moment  needs  to  be  numerically  calculated,  i.e.,  close 
form  solutions  such  as  that  for  the  box  girder  are  not  available. 

In  general,  the  fully  plastic  limit  state  is  not  useful  in  a  practical  sense  as  the  physical 
condition  it  represents  is  seldom  realized  because  of  buckling.  It  has  been  historically 
used,  however,  as  a  baseline  value  to  which  a  buckling  knockdown  factor  was  i^lied  in 
order  to  obtain  the  collapse  moment  for  the  hull  cross  section,  particularly  if  the  cross 
section  is  multicellular.  For  unicellular  cross  sections,  a  nxve  tq)propriate  baseline  value 
is  given  by  the  first  yield  moment  In  current  inactice,  the  buckling  knockdown  factOT  is 
applied  to  the  initial  yield  moment  as  indicated  in  Part  2  of  this  report 

10.1.3  Hull  Girder  Collapse  Limit  State 

The  first  yield  limit  state  and  the  fully  plastic  limit  state  are  both  idealizations  of  hull 
girder  behavior.  In  reality,  as  the  externally  applied  curvature  (or  moment)  on  the  hull 
girder  is  increased,  strains  internally  will  increase  up  to  a  point  where  either  the  yield 
strength  of  the  material  is  reached,  or  buckling  occurs  depending  on  the  slenderness  of 
the  structure.  Of  particular  importance  in  longitudinally  framed  vessels  is  the  buckling 
and  post  buckling  behavior  of  longitudinals  together  with  associated  plating,  and  also  in 
some  cases  the  overall  buckling  of  the  gross  panel  consisting  of  longitudinals  together 
with  the  transverse  beams.  When  parts  of  the  hull  buckle,  any  additional  load  is  "shed" 
to  or  taken  by  adjacent  stable  material,  up  to  the  point  at  which  they  also  buckle  or  reach 
yield.  As  the  externally  applied  curvature  increases,  typically  the  internal  resisting 
moment  calculated  with  accounting  of  buckling  and  yielding  in  parts  of  the  cross  section 
will  increase  up  to  a  point,  after  which  it  will  drop.  The  maximum  internal  resisting 
moment  so  calculated  is  the  so-called  collapse  moment,  'M^'.  On  the  tension  side  of  the 
hull  girder,  the  unloading/load  shedding  is  slower  and  on  the  compression  side,  it  is  more 
rapid.  A  typical  moment-curvature  diagram  for  a  hull  cross  section  is  illustrated  in 
Figure  10.3. 

We  have  not  speci^cally  considered  plate  buckling  in  the  above  discussion.  Buckling 
of  plate  elements  in  longitudinally  framed  situations  affects  the  collapse  moment  to  the 
extent  such  buckling  reduces  the  effective  width  of  plating  acting  with  longitudinals.  In 


91 


transversely  framed  situations,  the  plate  effect  on  collapse  mmnent  is  comparatively 
greater. 

CalciilatiQn  of  Collapse  Moment 

There  are  various  possible  methods  for  calculating  the  collapse  moment  These  vary 
from  iqtproaches  where  any  reserve  of  stiffened  plate  compresave  strength  after  its 
maTimum  resistance  has  been  reached  is  neglected,  to  ncnilinear  finite  element 
calculadmis  which  include  plastificatitxi  and  buckling  in  a  rigorous  way.  The  ctmeept  of 
downrating  or  knocking  down  the  fully  plastic  collapse  mranent  to  account  frv  buckling 
was  suggested  Caldwell  [IQ.  It  has  been  further  developed  by  Mansour  [7],  but  with 
knock  down  factors  to  be  qtplied  to  die  initial  yield  nxnnent  Procedures  incorporating 
an  incremental  moment-curvature  approach  to  hull  collapse  strength  have  been  developed 
by  Smith,  Billingsley  [17]  and  Adamchak  [18].  Finite  element  calculations  for  ship  hull 
colli^se  strength  are  presented  in  Thayamballi  et  al.  [19]. 

It  is  not  the  intendon  to  review  the  different  methodologies  for  ship  hull  collapse 
strength  calculation,  but  we  introduce  in  brief  here,  the  incremental  moment  curvature 
approach.  In  diis  method 

(i)  A  curvature  is  appUed  to  the  hull,  and  increased  incrementally. 

(ii)  Fev  each  value  of  curvature,  the  internal  resisting  moment  is  concqiuted, 
accounting  for  the  end  shemening  of  the  element  resulting  frtnn  internal 
strains,  including  any  buckling  and  post  buckling,  as  well  as  load  limitation  by 
plasticity.  Such  informatitm  is  included  through  load-end  shortending  curves, 
an  example  of  which  is  shown  in  Figure  10.4. 

(iii)  A  moment  curvature  relationship  for  the  hull,  such  as  that  in  Figure  10.3,  is 
developed,  and  the  collapse  moment  identified. 

The  most  important  part  of  the  calculations  is  the  establishment  of  the  load-end 
shortening  relationships  for  the  hull  members,  considering  the  various  local  failure 
modes. 
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lOJ  Limit  States  Associated  with  Local  Buckling 
10.2.1  General 


As  previously  noted,  these  define  the  strength  assnriatrri  with  column  and  beam 
column  buckling  of  longitudinals  together  with  associated  plating,  tripping  of 
longitudinals,  and  the  grillage  buckling  of  longitudinals  together  with  transverse  beams. 
The  strengths  calculated  do  not  account  for  any  post  buckling  reserve  which  is  typical 
small  (but  existent)  in  the  failure  modes  noted.  Also,  the  term  "local"  is  used  as  a 
qualifier  to  the  extent  that  only  one  component  is  considered  in  the  limit  state.  In  the  real 
structure,  there  may  be  several  such  identical  components  under  nominally  identical 
loading. 


10.2.2  Column  and  Beam  Column  Buckling 

Column  buckling  refers  to  the  flexural  buckling  of  longitudinals  together  with 
effecting  plating.  The  longitudinals  and  plating  may  be  pan  of  a  stiffened  panel  between 
transverse  beams.  The  panel,  and  hence  the  longimdinal  and  plating  are  considered  to  be 
under  compression.  In  the  beam-column  failure  mode,  in  addition  to  the  axial  load,  there 
are  also  lateral  loads  present  This  latter  situation  occurs  for  example  in  the  case  of 
longitudinals  and  plating  at  the  vessel  bottom.  The  column  idealization  is  shown  in 
Figure  10.5. 

Column  buckling  strength,  without  consideration  of  lateral  pressure,  is  given  by  the 
following  (Mansour,  Ref.  8): 


tt^E 

{Ijrf 


if  a„^ap 


o„>ap 


whCTe 


Cs 


ap(aY-ap) 
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The  first  equation,  valid  in  the  range  of  0  ^  ^  Cp,  where  Op  is  the  propoitional 

limit  stress  and  is  the  critical  buckling  stress,  will  be  recognized  as  the  Euler  elastic 
column  strength  equation.  In  the  second  equation,  a  correction  is  made,  based  on  a  factor 
Q,  if  the  calculated  elastic  buckling  stress  exceeds  the  proportional  limit  stress  Op.  The 
correction  is  such  that  the  limit  state  strength  calculated  from  the  pair  of  equations  given 
will  not  exceed  the  material  yield  strength. 

Also,  is  the  effective  column  length,  which  in  continuous  structures  where  the 
stiffener  ends  are  capable  of  rotation,  may  be  taken  equal  to  the  physical  length  between 
transverse  supports,  and  V  is  the  radius  of  gyration  of  the  cross  section  consisting  of 
plating  and  stiffener.  The  value  of  r  is  given  by: 


r  = 


where  1  tmd  A  tffe  the  mmnent  inertia  and  area  erf*  the  cross  seetkm,  re^)ectively. 
Typically,  in  computing  these  quantities,  an  effective  plate  flange  assuming  that  the  plate 
has  buckled  is  used.  The  plate  flange  width  may  be  obtained,  for  short  edge 
compression,  from  Mansour,  Ref.  8,  as  follows: 


b 


hi 

P 


if  PS3.5 


2.25  1.25 

P  P^ 


1.0  <p  <3.5 


=  1.0 


P<1 


where  P  is  the  non-dimensional  plate  slenderness,  defined  as: 


P  = 


b 

t 


where  b  is  the  width  of  the  short  edge,  i.e.,  the  spacing  of  longitudinals. 

In  the  case  of  beam-column  buckling,  the  lateral  pressure  results  in  a  reduction  in  the 
critical  stress  to  a  value  less  than  that  obtained  for  the  column  buckling  limit  state.  A 
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relatively  siiiiple  apjffoach  to  characterizing  limit  state  strength  for  this  situadcm  is  to  use 
a  linear  interaction  equation: 


— J-+— i.  =  l 
Oc  Oy 

where  ^  ^  arfumn  buckluig  stteagtb  assummg  no  btten^  pressure,  aai  is  the 
yield  strength,  a,  is  the  axial  stress  and  Oj,  is  the  maximum  bending  stress  over  the  span 
of  the  ItmgitudinaL  This  interaction  equation  assumes  duu  tripping  of  the  cross  sectim, 
and  local  buckling  in  the  cross  section  (e.g..  of  the  flange  or  web)  are  avoided. 

The  calculation  of  o,  should  account  for  any  reductkm  in  plate  effectiveness  because 
of  buckling.  The  calculation  of  the  bending  stress  should  in  principle  account  for  shear 
lag  effects,  although  for  panels  with  closely  spaced  Imgitudinals,  the  effect  may  often  be 
neglected. 

10.2.3  Tripping  of  Longitudinals 

In  this  failure  noode,  also  called  torsional/flexural  buckling,  failure  is  initiated  by 
twisting  of  die  stiffener  in  such  a  way  dot  dn:  joint  between  die  stiffener  ami  phtte  does 
not  move  laterally.  A  portion  of  the  adjacent  plate  may  participate  in  the  twisting,  and 
the  flange  of  the  stiffener  may  twist  together  with  the  web,  or  the  two  may  twist 
differentially.  Tripping  is  illustrated  in  Hgure  10.5.  The  trij^ing  phenomenon  may 
occur  under  axial  loads  alone,  or  under  axial  loads  in  combination  with  lateral  pressure 
loads. 

The  ultimate  strength  for  torsional/flexural  buckling  under  axial  conqiressive  loading 
may  be  obtained  as  follows  (see  Reference  8): 

a)  Calculate  the  elastic  tripping  stress  a,  for  the  stiffener  cross  section  rotating  about 
an  enforced  axis  at  its  toe.  This  is  given  by 


where  G  is  the  shear  tnodulus  for  the  material,!  is  the  tonaott  constant,  and  is 
the  warping  constant.  The  length  of  the  longitudinal  between  supports  is  denoted 
T.  Expressions  for  the  torsion  and  warping  constants  as  a  function  of  cross 
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section  shi^  may  be  found  in  the  book  by  Bleich,  Ref.  20.  is  the  polar 
.  moment  of  inertia  of  the  cross  section  about  an  enfoiced  axis  at  its  toe.  Le., 

Io»I,+I,  +  Aya 

where  and  I,  are  the  principal  moments  of  inertia  of  the  cross  section,  of  area 
A,  and  y  is  the  web  depth. 

b)  Obtain  the  elastic  tripping  stress  considering  interaction  with  column 
buckling,  by  solving  the  following  quadratic: 

o.  *  0 

*o 

Here,  1^  is  the  polar  moment  of  inertia  of  the  cross  secdM  of  the  stiffener,  i.e., 

+  ly,  and  ocr  is  the  limit  state  strength  fen-  colunm  buckling  under  axial  loads.  If 
otfe  ^  Op,  where  <Sp  is  the  proportitmal  limit  stress,  the  tripping  limit  stress  « 
Otherwise,  otf  is  obtained  from 


The  above  determination  of  limit  state  strength  for  tripping  of  longitudinals  under 
axial  loads  is  outlined  in  Mansour,  Ref.  8.  When  lateral  pressure  is  present,  the  axial 
tripping  strength  should  be  modified  to  reflea  its  influence.  Although  more  detailed 
approaches  are  possible,  one  way  to  inclixle  the  lateral  pressure  ejects  is  to  use  a  linear 
interaction  formula  similar  to  that  used  for  the  case  of  the  beam-column  limit  state.  Such 
an  approach  will  not  apply  to  a  case  where  the  pressure  loads  are  the  dominant  ones,  and 
additional  refinements  will  be  needed. 

10.2.4  Grillage  Buckling 

This  failure  mode  and  the  limit  state  strength  associated  with  it  refer  to  the  buckling 
of  the  gross  panel,  i.e.,  longitudinals  and  transverses,  between  the  maja*  suppon 
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members  such  as  bulkheads.  A  portion  of  such  a  gross  panel  under  compression  is 
shown  in  Figure  lO.S.  This  problem  has  been  extensively  studied  by  Mansour  [21,22] 
using  orthotropic  plate  theory.  The  following,  taken  from  Ref.  22,  may  be  used  if  the 
number  of  stiffeners  in  each  direction  is  sufficiently  large,  e.g.,  3  to  S. 

For  gross  panels  under  uniaxial  compression,  the  critical  buckling  stress  is  given  from 

“  ■  h.  B’ 

where  B  is  the  width  of  the  gross  panel,  h,  is  the  effective  thickness  resisting  the 
compressive  loads  in  the  x  direction,  and  k  is  a  buckling  coefficient  that  depends  on  the 
boundary  conditions.  For  simply  supported  gross  panels. 


k  = 


m*  ^ 
p  m 


For  gross  panels  with  both  loaded  edges  simply  supponed  and  both  the  other  edges  fixed. 


2  *2 

k  =  ^+2.5Tl  +  5-iV 

p‘‘  m 


where  m  is  the  number  of  half  waves  of  the  buckled  orthotropic  plate,  to  be  chosen  such 
that  k  is  minimized;  T)  and  p  are  the  virtual  aspect  ratio  and  the  torsion  coefficient, 
respectively. 

The  virtual  aspea  ratio  and  the  torsioa  coefficient  are  given  by 


n= 


P  = 


L 

B 


Here  (see  Figure  10.6),  and  D  are  the  flexural  rigidities  per  unit  width,  given  by 


D  = 


El. 


’  S,(l-v=) 


;  D,= 


El. 


S.(l-v^) 
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where  I,  and  1,  are  the  nKxnents  of  inertia  of  the  stiffeners  extending  in  the  x  and  y 
direcdtMis  (i.e.,  about  the  y,  x  axes),  and  Ip^  are  the  moments  of  inertia  of  the 
effective  plate  flange  alone,  acting  with  the  stiffeners  in  the  x,  y  directions.  and  Sy  are 
the  X  and  y  stiffener  spacings. 

The  effective  plate  thickness  h,  is  the  average  cross  sectional  area  per  unit  width  of 
effective  plating  and  stiffeners  in  the  x  direction,  i.e. 

.  A, +3,  t 

’’■'“s — 

y 

where  A,  is  the  stiffener  area,  t  is  the  plate  thickness,  and  is  the  effective  width  of  the 
plate  flange,  ^  Sy. 

Reference  21  by  Mansour  contains  an  extensive  treatment  of  the  behavior  of 
orthotropic  plate  panels  in  the  buckling  and  elastic  post  buckling  range.  Design  charts 
are  given,  which  address,  for  example,  the  midplane  deflection,  critical  buckling  stress, 
and  the  bending  moment  at  midlength  of  the  edge.  The  types  of  loading  considered 
include  combinadons  of  normal  pressure,  direct  inplane  stresses  in  two  directions,  and 
edge  shear  stress.  From  the  charts,  prediction  of  large  deflection  behavior  up  to  the  onset 
of  yielding  is  possible  in  a  practical  sense.  Alternatively,  in  a  unidirectional  load 
situation  which  is  a  very  common  case,  Ihnit  state  strength  may  be  obtained  from  the 
previously  given  close  form  expressions  from  Ref.  22.  That  solution  is  not  valid  beyond 
the  linear  elastic  regime,  unless  corrections  of  the  type  made  for  column  behavior  are 
also  made  in  this  case. 
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Figure  10.2  Development  of  the  Fully  Plastic  Limit  State 


Figure  10.3  Moment-Girvature  Diagrams  for  a  Ship  Hull 


Figure  10.4  Load-End  Shonening  Curve  for  a  Column 
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Figure  10.5  Stiffener  Plate  Failure  Modes 
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11.  Conclusions  and  Discussion 


11.1  Summary  and  Major  Results 

Two  demonstrations  have  been  carried  out  in  this  project;  a  demonstration  of 
probability-based  Rule  calibration  (Part  1),  and  a  demonstration  of  probability-based  hull 
girder  safety  analysis  (Part  2).  Also,  an  extension  to  the  project.  Part  3  defined  loads, 
strength  and  structural  reliability  terminology,  identified  ultimate  and  serviceability  limit 
states,  and  considered  procedures  for  load  extrapolation  and  load  definition. 

In  the  first  part,  the  calibration  procedure  was  described  and  applied  to  ABS  hull 
girder  longitudinal  strength  formulation.  For  this  purpose  300  "ABS  Ships"  are 
considered  and  the  minimum  required  section  modulus  of  each  has  been  determined 
according  to  ABS  Rules  (see  Appendices  1  and  2).  The  safety  index  p  was  then 
determined  using  first  and  second  order  reliability  methods.  It  was  found  that  the  safety 
indices  vary  slightly  and  that  variation  depended  only  on  the  ratio  of  the  wave  bending 
moment  to  the  Stillwater  bending  moment.  The  range  of  the  safety  indice.s,  Prangc'Pnuix" 
Pniin*  found  to  be  0.31.  The  average  value  of  the  safety  indices  was  found  to  be 
3.2. 

The  aim  of  the  calibration  procedure,  which  is  described  in  detail  in  Part  1  of  the 
report,  is  to  eliminate  this  variation  in  §  in  order  to  achieve  uniform  safety  standard  for 
all  ship  sizes.  The  target  p  value  was  taken  as  the  average  value,P(=fl,y=:3.2.  The 
calibrated  formulation,  which  is  based  on  partial  safety  factor  format,  produced  the  target 
value  of  p  and  a  Pj,„ge=0.004, 

It  should  be  noted  that  the  calibrated  formulation,  in  as  much  as  the  initial  ABS 
formulation,  ensures  :)ni  ■  a  safely  level  against  deck  yielding.  For  buckling 
considerations,  the  stiffening  system  for  each  of  the  300  "ABS  Ships"  must  be  designed 
and  evaluated.  Buckling  rule  calibration  is  best  done  at  the  local  level  since  the  Rules 
control  and  specify  stiffener  spacing,  section  modulus  and  plate  thickness  at  a  local  level. 
Similar  calibration  procedure  to  that  described  in  Part  1  can  be  used  to  calibrate  ABS 
formulations  that  give  minimum  required  stiffener  section  modulus  and  plate  thickness  so 
as  to  produce  uniform  safety. 

In  Part  2  of  the  report  a  tanker  was  token  as  an  example  to  demonstrate  the  use  of 
probability-based  safety  analysis,  i. '  ef'timate  the  reliability  in  an  existing  ship  (or  on 
a  drawing  board  design).  For  this  purpose  several  limit  states  have  been  developed 
including  ultimate  strength  (buckling  collapse,  deck  initial  yield  and  fully  plastic 
collapse),  serviceability  limit  state  (local  plate  buckling)  and  fatigue  limit  state.  Mote 
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realistic  load  estimates  have  been  developed  for  each  limit  state,  based  on  parametric 
seakeeping  and  ship  motion  analysis.  The  wave  bending  moment  has  been  calculated  for 
the  ultimate  limit  state  with  considerations  given  to  the  most  probable  extreme  sea 
condition  the  ship  is  likely  to  encounter.  For  the  fatigue  limit  state,  stress  ranges  and 
number  of  cycles  have  been  calculated  based  on  a  sea  scatter  diagram. 

A  reliability  index  p  has  been  calculated  using  first  and  second  order  reliability 
methods  for  each  limit  state.  Model  uncertainty  was  included  in  all  limit  states.  The 
resulting  safety  indices  indicate  that  buckling  collapse  is  the  governing  mode  of  failure  as 
its  safety  index  is  well  below  those  of  deck  initial  yield  and  fully  plastic  collapse. 


11.2  Benefits  and  Drawbacks  of  Using  Probability -based  Design  Method 

Use  of  probabilistic  methods  in  design  can  provide  several  benefits  and  some  unique 
features.  Among  those  are: 

1.  Explicit  consideration  and  evaluation  of  uncertainties  associated  with  the  design 
variables. 

2.  Inclusion  of  all  available  relevant  information  in  the  design  proce.ss. 

3.  Provides  a  framework  of  sensitivity  measures. 

4.  Provides  means  for  decomposition  of  global  safety  of  a  structure  into  partial  safety 
factors  associated  with  the  individual  design  variables. 

5.  Provides  means  for  achieving  uniformity  of  .safety  within  a  given  class  of  structures 
( or  specified  nonuniformity ). 

6.  Minimum  ambiguity  when  updating  design  criteria. 

7.  Provides  means  to  weigh  variables  in  terms  of  their  significance. 

8.  Provides  rationale  for  data  gathering. 

9.  Provides  guidance  in  novel  design. 

10.  Provide  the  potential  to  reduce  weight  without  loss  of  reliability,  or  improve 
reliability  without  increasing  weight.  The  methods  can  identify  and  correct  overly 
and  unduly  conservative  designs. 


106 


In  addition  to  the  above  benefits,  reliability  technology  lends  itself  for  certain  use  for 

which  it  is  much  more  suitable  than  traditional  design  methods.  In  reference[14], 

Wirsching  lists  some  of  its  use,  which  include: 

1.  To  compare  alternative  designs,  particularly  in  the  early  stages  when  several 
competing  design  concepts  are  considered. 

2.  To  perform  failure  analysis  of  a  component  or  a  system. 

3.  To  develop  a  strategy  for  design  and  maintenance  of  structures  which  age  (e.g., 
corrosion,  fatigue),  and  to  determine  inspection  intervals. 

4.  To  execute  "economic  value  analysis"  or  "risk  based  economics"  to  produce  a  design 
with  a  minimum  life  cycle  costs. 

5.  To  develop  a  strategy  for  design,  warranties,  spare  parts  requirements. 

6.  In  general,  as  a  design  tool  to  manage  uncertainty  in  engineering  problems. 


Use  and  implementation  of  probabilistic  methods  are  not  without  problems.  Some  of 
the  drawbacks  are: 

1.  Use  of  reliability  analysis  in  safety  and  design  processes  requires  more  information  on 
the  environment,  loads  and  the  properties  and  characteristics  of  the  structure  than  typical 
deterministic  analyses.  Often  .some  information  are  not  available  or  may  require 
considerable  time  and  effort  to  collect.  Time  and  schedule  restrictions  on  design  are 
usually  limiting  factors  on  the  u.se  of  such  methods. 

2.  Application  of  probabilistic  and  reliability  methods  usually  require  some  familiarity  of 
basic  concepts  in  probability,  reliability  and  statistics.  Practitioners  and  designers  are 
gaining  such  familiarity  through  seminars,  symposia  and  special  courses.  Educational 
institutions  are  also  requiring  more  probability  and  statistics  courses  to  be  taken  by 
students  at  the  graduate  and  undergraduate  levels.  This,  however,  is  a  slow  process  that 
will  take  some  time  in  order  to  produce  the  necessary  "infrastructure"  for  a  routine  use  of 
reliability  and  probabilistic  methods  in  design. 

3.  On  a  more  technical  aspect,  the  reliability  analysis  did  not  deliver  what  it  initially 
promised,  that  is,  a  true  mea.sure  of  the  reliability  of  a  structure  by  a  "true  and  actual" 
probability  of  failure.  Instead  what  it  delivered  is  "notional  probabilities"  of  failure  and 
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safety  indices  which  are  good  only  as  comparative  measures.  Only  notional  values  are 
delivered  because  of  the  many  assumptions  and  approximations  made  in  the  analysis 
producing  such  probabilities  and  indices.  These  approximations,  deficiencies  and 
assumptions,  however  are  made,  not  only  in  probability-based  design,  but  also  in 
traditional  design.  Approximations  are  made  in  the  determination  of  loads  using 
hydrodynamics  theory  and  in  the  structural  analysis  and  response  to  the  applied  loads. 
When  all  such  assumptions  and  deficiencies  are  removed  from  the  design  analysis,  the 
resulting  probabilities  of  failure  will  approach  the  "true"  probabilities. 


lU  Discussion  of  SSC  Projects  in  Reliability  and  Needs  to  be  addressed  in 
Further  Projects 

The  strategic  plan  of  the  Committee  on  Marine  Structures  (CMS)  as  outlined  in  the 
Marine  Board  report  entitled  "Marine  Structures— Research  Recommendations  for  FY 
1992"  has  been  reviewed.  In  this  document,  the  CMS  states  the  goals  and  objectives  of 
the  plan  and  lists  a  five-year  research  program  and  development  which  is  organized 
under  five  technology  areas.  The  technology  areas  are:  reliability,  loads  and  response, 
material  criteria,  fabrication  and  maintenance,  and  design  methods.  The  five  technology 
areas  consist  of  23  compreheasive  and  well  thought-out  subject  areas.  The  projects 
outlined  in  these  subject  areas  will  undoubtedly  lead  towards  fulfilling  the  goals  of  the 
plan  which  include  improving  the  safety  and  integrity  of  marine  structures,  improving 
competitivene.ss  of  U.S.  merchant  shipping,  and  promoting  the  development  of  new 
marine  systems. 

Based  on  the  work  carried  out  in  this  project  and  the  review  of  CMS  research 
recommendations,  the  following  areas  are  suggested  for  further  development.  Some  of 
these  areas  are  very  specific  and  each  need  to  be  addres.sed  in  depth  as  a  limited  scope 
project.  These  gaps  are: 

1.  Torsional/flexural  buckling  (tripping  failure)  of  ship  stiffeneis  with  effective  breadth 
of  plating  —  analysis  and  development  of  design  formulation. 

2.  Ultimate  strength  of  ship  hull  girders  due  to  instability  -  analysis  to  determine 
strength  reduction  factors  due  to  instability  to  be  used  in  design. 
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3.  Experiments  on  hull  girder  ultimate  strength  to  verify  analytically  calculated  strength 
reduction  factors. 

4.  Selection  of  wave  spectra  (or  wave  data)  pertinent  to  design  wave  loading  on  ships. 

5.  A  study  leading  to  the  determination  of  the  ratio  sag  to  hog  wave  bending  moments 
and  the  bias  associated  with  linear  ship  motion  load  prediction. 

6.  Design  formulation  for  combined  wave  and  slamming  bending  moments. 

7.  A  study  of  shear  forces  and  moments  acting  on  the  forward  part  of  a  ship  including 
slamming  effects. 

8.  A  study  leading  to  target  reliabilities  for  each  hull  girder  limit  state  ba.sed  on  existing 
ships. 

9.  Development  of  reliability  procedures  and  target  indices  for  local  structure  in  ships. 

10.  A  study  to  develop  a  re  I  lability -batted  cost  analysis  which  aims  to  achieve  minimum 
life  cycle  costs  for  ships. 

11.  Development  of  a  reliability-bH.sed  strategy  for  inspection  intervals  and  maintenance 
of  ships. 

12.  Inclusion  of  system  reliability  considerations  in  fatigue  and  multiple  failure  modes. 

13.  Reliability  assessment  of  transverse  structures  and  lateral  pressure  effects. 
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APPENDIX  1 

Msw,  Mw,  Mw/Msw,  and  SM  of  "ABS  Ships" 
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APPENDIX  2 

Means  and  Standard  Deviations  of  Msw,  Mw,  and  SM  of  "ABS 
Ships" 
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APPENDIX  3 

3.1  Calculation  of  Plastic  Moment  Capacity 

3.2  Calculation  of  Critical  Buckling  Stresses 

3.3  Calculation  of  Effective  Section  Modulus  after  Buckling 
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3.1  FULLY  TLASTIC  MOMENT  CAPACITY 


Mp  =  fully  plastic  mcHnent  =  (SMp)  ■  fy 

fy  =  yield  strength  of  the  mateiial  =  259  N/rnm^  (37.6  ksi) 

(SM)p  =  plastic  section  modulus 

From  SSC  219  "Ultimate  Strength  of  a  Ship's  Hull  Girder  in  Plastic  and  Buckling 
Modes": 

(SM)p  =  Ajjg  +  2(A5  +  Ablk)  ^  ■  8  +  ■  8) 

J.  Afl  +  2(As  +  Ablk)  .  Ap 

D  ^  4As  “ 

D  =  24  m  =>  g*  14.181m. 


(SM)p  =  5.8376  •  10^  mcm^ 


Ratio  between  plastic  section  modulus  and  the  elastic  section  noodulus: 


(SM)p  5.8376  •  IQS  _ 

(SM)c  ^  4.65767  •  KP  " 

Also, 

Aj)  =  1.4645  •  10®  mm2 

Ab  =  1.9934- 10®  mm2 

As  =  7.9654  •  105  mm2 

Ablk  =  6.5830  •  105  mm2 
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TANKER  MIDSHIP  SECTION 

Scantlings  called  out 
where  used  in  calculations 


Calculations  follow  Ref.  [8] 


Considering  only  verdcal  bending  moment>  so  uniaxial  compressive  stress: 


^  2.25  1.25 


1.0 


if  ^  :^3.5 

I.O  <  P  <  3.5 

1.0 
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ServiceabUitv  Limit  State 


334.82 

155.4(259-155.4) 


6.96 


334.8 


iL 

nm^ 


>  a. 


155.4 


N 

nm* 


«J6(259)^^^_N 
•  6.96+1 _ TDB^ 


n.  Stiffeners  and  Effective  Plating 

The  compressive  strength  of  the  stiffeners  together  with  effective  plating  is 
considered.  Only  ultimate  limit  state  is  considered,  because  when  a  column  buckles  it 
reaches  immediately’ its  uhnnate  strength. 

The  effective  plating  under  edge  compression  is  determined  from; 


be 


=  1000  0.92  =  920mm 
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be2  (from  buckling  considerations)  =  0.221(5400)  =  1193  mm 
bc2  (firom  shear  lag  analysis)  =  0.9(5400)  =  4860  mm 


A 

Ix 

ly 

e 

lo 

Ic 


=  3.2816  •  10*  nim2 
=  3.2816  •  10*  mm* 

=  1.363  - 10®  mm* 

=  450.5  distance  from  neutral  cxis  to  shear  center  =  96.9  mm 

=  lx  +  Iy  +  V  =  2.34  •  10®  mm* 

=  I,  ly  =  2.035  •  10®  mm* 


J  s  torsional  const 


920  -  213 +  (450+^1  303 

= - 4 ^ - =  7.08  • 


10®  mm* 


21-4502  9203-303  . . 

=  warpmg  constant  =  9203  +  303  ® 


2.1  - 103 


N 


^=i(GI.’^)=  244£5 


i)  Elastic  Range:  Consider  interaction  with  flexural  buckling. 

I  N 

f  -  <3tfe'<3tfe(Ocr  +  <3t>  +  <3er-<3t  =  0  nim2 

*a 
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ii)  Plastic  Range: 


( 


1  - 


Uniaxial  Compressive  Load  -  Serviceability  Limit  State 
C,  =  2662.6  mm 


I,  =  1.4386  •  IQio  +  4.44  •  lO^  +  1.133  •  10“  =  1.3213  •  10“  mm^ 

Ip,  =  1.4386- 10»>mm« 

For  the  calculation  of  ly  and  Ipy  an  effective  breadth  of  4860  is  used: 

Cy  =  1357.33  mm 

ly  =  4.072  •  109  +  3.25  -  lO^o  =  3.6572  •  lO^o  irnn^ 

Ipy  =  4.072- lO^mm^ 

Sy  =  1000  mm 
S,  =  5400  mm 

t,  =  equivalent  thickness  of  plate  and  stiffeners  extending  in  x-direction 


450-30-18 2850-18 
20000 


•J-21  =  35.715  mm 


n 

"  S„(l  -  v2) 


=  3.049  - 10*3 


D„  = 


Ely 


y  S,(l  -  v2) 


=  1.563  - 10^2 


A/B  =  1.08 


V  =  0.3 
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9(6.692)  •  10*  =  6.0228  •  lO’ 


ly  =  3.6572- 10^0  mm* 
Sy  =  1000  mm 
Sx  =  5400  mm 


18 


From  buckling  considerations  =0.0597  -  20  =  1193  mm 
A  =  7.635  - 10*  mm2 
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C  s  1896.1  mm 


=  6.94  •  10*®  mm* 


nf  ai  ai  XU  M  Vi  Dici  d[m  i  ciici/i  (ii  uii  nik'^i  ai  m  ai ; 


inicLMicioxM; 


b  =  92%  of  original  width 


C  s  distance  from  local  neutral  axis  to  global  neutral  axis 
I  =  4.657675  •  10*®  mm  •  12950  =  6.0137  •  10*^  mnri 

SMeff  =  j2950 ^ ^  =  4.570443  •  10*®  mm^ 


|SH.ff  =  4.570443  •  lO^mcm^l  reduced  1.9% 
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APPENDIX  4 

Calculations  of  Compressive  Strength  Factor  and  the  Hull  Girder  Instability 
Collapse  Moment 
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The  CoQ^iressive  Strength  Factor  fw  the  Critical  Panel  of  the  Example  Ship 
(ISSC  Formula) 


9q) 

X 

P 


=  (0.960  +  0.765  X.2  +  0.176  +  0.131  X^p  +  1.064  X^)^^ 

_  X.  £  _  ^400-- inr 

“  m\E  ~  142.8  irV 2.1  • 


105 


=  0.403 


£  1000-  /  235Z 

*7  VE  =  21  V2.1.105 


=  1.593 


=  0.787 

For  the  sagging  condition,  we  then  have 

Mu  =  (-0. 172  +  1 .548  9q,  *  0.368  9^  SM  •  fy 
=  0.819  SM- 4 
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APPENDIX  5 

Calculations  of  the  RMS  Values  of  the  Wave  Bending  Moment  for  the 
Example  Ship 

5.1  Ultimate  Limit  State 

5.2  Fatigue  Limit  State 
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5.1  RMS  OF  EXTREME  WAVE  BENDING  MOME^ 
(ULTIMATE  LIMIT  STATE) 

VESSEL  AND  SEA  STATE  DATA 
Cb  »  0.71 
L/B  =  6.19 
B/T  =  2.62 

CALCmJlTlONPROCEIXJRE 

Calculations  are  according  to  seakeq)ing  tables  of  Ref.  6.  From  the 
seakeeping  table  (see  sample  interpolation  chart  on  the  next  page), 

rms  =  272.7 

This  value  is  made  dimensional  by  multiplying  it  with:  pg 

where  p  «  specific  density  of  seawater  «  1025  kgAn^ 
g  =  acceleration  of  gravity  =  9.81  m/s^ 

L  s  length  of  ship  »  260  m. 

Dimensional  tins  s  125398'  lO^kNm 

This  value  may  be  overestunated  a  few  percent  because  a  Roude  number  of  0.1  is 
s^lied  instead  of  the  value  0.05. 

The  seakeeping  tables  are  not  tabulated  for  values  of  lower  than  0.1. 


Hs  =  12.2  m  (40  ft) 

S  =  Hs/L  =  0.047 

F„  =  0.05  (will  use  F„  =  0.1) 
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INTCMOtATION 
TABLE  II 


£  * 
N  s 


w  . 

Xa  I  1 

-  «n 
21 


o 

rt> 

m 


^i-f' 

Lcir'“ 

-‘IL. 

=Jr” 


cnI _  ♦ 

■=  Jr 

^1L.. 


♦_  s 


-  >s 


'-1 

zF-  IL  u 


T' 

Imwm» 


N, _ ♦ 

z  Jr 

^iL. 


J  o 


C*" 


•”-c: 

♦ 

f..M  ■ 


hi: 


Interpc^tion  chart 

S— k—ping  Standard  Series  for  Crul8er*Stern  Ships 


S.2  RMS  VALUE  FOR  WAVE  BENDING  MOMENT  (FATIGUE  UMIT  STATE! 


Hs[m] 

rms  [kNm] 

0.5 

3.1705  •  10* 

1.5 

9.6541  •  10« 

2.5 

1.6639  •  105 

3.5 

2.1385  •  105 

4.5 

3.3420  •  105 

5.5 

4.8565  •  105 

6.5 

6.2111-  105 

7.5 

7.4853  •  105 

8.5 

7.9416  •  105 

9.5 

9.5985  •  105 

10.5 

1.0340  •  10^ 

11.5 

1.1082-105 

12.5 

1.1686-105 

13.5 

1.2404  -  105 

The  above  results  are  for  the  sea  scatier  diagram  used  in  the  fatigue  analysis  and 
shown  in  Appendix  F.  The  inteipolatioD  charts  using  the  seakeeping  tables  of  Ref.  6  are 
omitted  for  brevity,  but  each  calculation  is  amilar  to  that  previously  shown  for  the  rms  of 
extreme  wave  bending  moment 
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APPENDIX  6 

Fatigue  Reliability  CakubtioQS 

6. 1  Fatigue  Reliability  Analysis  of  I>eck  Detail 

6.2  Sea  Scatter  Diagcama 
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6.1  FATIGUE  RELIABILITY  ANALYSIS  OF  DECK  DETAIL 


The  detail  is  shown  in  Figure  6.1  and  classified  as  belonging  to  class  D  [13].  The 
long  tenn  statistics  of  sea  states  is  from  the  Oseberg  Area  of  the  North  Sea.  It  is  shown 
elsewhere  in  this  section. 


The  class  D  gives  the  S-N  curves: 

log  N  =  log  a  -  2  log  s  -  m  log  AS 

=  1 1.7525  -  2  ■  0.1793  -  3  •  log  AS 

N  =  number  of  cycles 

A(S  =  stress  range 


C 


=  NAS®  =  100^6007 -2  0.4190)  =  1.52  •  10*2  N/inm2 


The  limit  state  function  is 


^ .  c 

T  is  the  service  life  of  the  ship  =  20  years. 

Q  is  the  stress  parameter  which  is  given  below: 


D  = 


{2-J2T 

2k 


(1) 


where 

m  =  fixed  =  3  (from  SN-curve) 


Xgj,  X2j  are  zero  and  second  stress  spectrum  moment  in  j-th  sea  state. 

From  the  seakeeping  tables  [6],  the  nns  for  the  wave  bending  moment  is 
obtained.  The  relation  between  the  zero  stress  spectrum  moment  and  zero  wave  b.m 
spectrum  moment  is: 


6-2 


^stnsis  f  ^  distance  from  NA  to  fatigue  cntckT  ^  wa^ 
o  “  distance  from  NA  to  deck  J  ° 

For  the  example  ship; 

>  itrw  =  (4.2948  •  1(H  [m6]-i>5f®M 


and 


5)Pj  •  X'jf  =  2-009-10^^[kNmP[secr^ 


when  theP^j  andX2j  ate  for  the  wave  bending  moment  obtained  in  Appendix  E. 
Equations  (1),  (2)  and  (3)  give  Cl,  the  stress  parameter; 


Cl  = 


(2\l2) 

2ir 


r  (4.2948  •  l(H)3/2  •  2.009  •  lO^^  [sec]-i 


(2) 


(3) 


=»  fl  =  852  [sec]*i 
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APPENDIX  7 

Typical  Input/Output  File  of  CALREL 

7.1  User  Defined  Subroutine  for  Limit  State  Function  and  Wave 
Bending  Moment  Distribution 

7.2  Input  Data  Hie 

7.3  Output  Hie 
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CALMl  nrxaS  ncp«l 
DATA 

TXTL  nlilM  CltX* 

1 

•xaapla  ship  rallabillty  Analysis  ••  dsck  initial  yiald.  casal 
PLMS  icl.igr 
1  0 

OPTZ  i^.nll.ni2.tol.opl.op2.op3 
1.20.4.0.001 

STAT  igt  (i)  .nga.nga  nv.ids,ax.sg.p3,p4.x0 

1  8 

SB  1.2.4.57a5.1.828a4 
Cp  2.2,25.8.1.813 
sw  3,1.3.022a6.1.0 

w  4.-51.4.855a6.4.3695aS.0.0,0.0,4.8S5a< 

XU  5,1.1.0.0.15 

xsw  6,1.1.0.0.05 

Xkr  7.1.0.9,0.135 

xs  8.1.1.15.0.0345 

END 

FORM 

SENS 

SORM 

EXIT 
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lapllcit  rMl*8  («-h.o*z) 
dlMIMiOR  x(l).tp(l) 

9  m  X(5)*X(1)*X(2)-X(C)*X(3)-X(7)*X(S)*X(4) 

raturn 

•nd 


•ubroucln*  udgx(dgx,x,tp.ig) 
laplicit  rMil«8  («-h,o-s) 
dlBMnslon  x(l),dgx(l).tp(l) 
return 
•nd 


•utnroutina  udd(x.p«r,sg,ids.ed£.pdf ,bnd,lb) 
iaplleit  r««l«8  (a-h,o-z) 
dlMnaion  x(l)  ,par(4)  .bnd(2) 
pi>3. 1415926 

fact 1> ( sqrt ( 6 ) /pi ) •par (1 ) *p«r (2 ) • <0 . 5772 *6/pi **2 ) * (par  (2 ) **2 ) 

£ace2 «d«xp ( 0 . 5 * ( (pi /aqrt ( 6 ) ) • (par ( 1 ) /par ( 2 ) ) -0 . 5772 ) ) 
cd£«daxp(-£act2*daxp(- (x(4) **2) / (2*£actl) ) ) 
pdf. (x(4) /factl) *fact2*daxp(-  (x(4) **2) / (2*factl) ) *cdf 
bnd(l).0.0d0 

Ibal 

«g>par(2) 

raeum 

and 
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7- 


Univarslty  of  C«llCernl« 
OttparCBMit  of  Civil  EngiiiMriag 

C  A  L  R  I  L 
CAL'RCLiability  prograa 
0«vttlap«d  by 

P.-L.  Liu.  N.-Z.  Lin  «nd  A.  Onr  Kiuraghian 

Lane  Ravition:  January  1990 
Copyright  0  1990 


MARNZNS  2<  coonand  not  availabla 


»»  NEW  PROBLZM  «« 


nunbar  of  liait-atata  ftinctiMu . ngf«  1 

nuabar  of  indapandant  variabla  groi^  ...nig*  1 

total  nunbar  of  randon  variablaa  . nrx*  • 

nunbar  of  linit-stata  paranatars  . ntpa  1 


»»  INPUT  DATA  «« 


axanpla  ship  raliability  analysis  —  dack  initial  yiald,  casal 

typa  of  systan . Ida  1 

iclsl  . conponant 

lels2  . sarias  systsn 

icla3  . .ganaral  systan 

flag  for  gradiant  conputation . igra  0 

igraO  . finita  diffaranea 

igral  . fomulas  providad  by  usar 


optinization  sehana  usad  . iopa  1 

iopai  . HL'RF  nathod 

lopa2  . nodifiad  HL^KF  nathod 

iopai  . gradiant  projaetion  nathod 

iopa4  . saquantial  quadratic  nathod 

naxlmun  nunbar  of  itaration  eyclas . nila  20 

HHUciimm  staps  in  lina  saareh . ni2a  4 

convarganca  tolaranca  . tola  l.OOOE-03 

optinization  paranatar  1  . opla  l.OOOE^OO 

optinization  paranatar  2  . op2a  O.OOOE^OO 

optinization  paranatar  3  . op3a  O.OOOEf^OO 


statistical  data  of  basic  variblas: 
availabla  probability  distributions: 


dataminitic  . ids«0 

nomal  . idssl 

lognonsal  . ids«2 

gasna  . ids«3 

shiftad  axponantial  . idss4 

shiftad  raylaigh  . idssS 

unifom . idssO 

bata  . ids«7 

typa  i  largast  valua . idssll 

typa  i  smallast  valua  — .ldsal2 
typa  ii  largast  valua  ....idssll 

waibull  . Idsal4 

usar  dafinad . ids>50 


group 

no. 

:  1 

var 

ids 

naan 

sn 

2 

4.S7E-»05 

Cp 

2 

2.59E401 

sw 

1 

3.02E406 

mu 

51 

XU 

1 

l.OOE-i'OO 

xsw 

1 

l.OOE-fOO 

xw 

1 

9.00E-01 

xs 

1 

1.15E4^00 

group  typa: 


St.  dav. 

paranl 

1.83E>04 

1.30E4-01 

1.81E4-00 

3.25E4-00 

1.00E<»00 

3.02E^06 

4.37E-»05 

4.86E406 

l.SOE-01 

l.OOE^OO 

5.00E-02 

l.OOE^OO 

1.35E-01 

9.00E-01 

3.45E-02 

l.lSE-»00 

1 


paran2 

paranl 

4.00E-02 

6.99E-02 

1.00E4-00 

4.37E+05 

0.00E*00 

l.SOE-01 

5.00E-02 

1.35E-01 

3.45E-02 

paran4 


O.OOE-fOO 


init.  pt 
O.OOE'fOO 
0.00E4>00 
O.OOE-fOO 
4.86E-f06 
0.0OB400 
O.OOE-i'OO 
O.OOE-fOO 
0.00E<<>00 


print  interval  . npr*  o 

npr<0  . np  first  erdsr  rssults  srs  priatsd 

npr-O  . print  ths  final  ttsp  of  FORM  rssults 

npr>0  . print  ths  rssults  of  svsry  npr  stops 

initial izst ion  flag  . inl«  0  7-5 

inisO  . start  froa  saan  point 

inial  . start  frosi  point  spscifisd  by  ussr 

inia-i  ....start  froa  prsvious  linsarisation  point 

rsstart  flag  . ista  0 

istaO  . analyzs  a  now  problsa 

istai  . continus  an  unconvsrgsd  problsa 


limit-stats  function  1 


itsration  nuabsr  . itsra  7 

valus  of  limit-stats  function. .g{x)a-2.805E-05 

reliability  index . bstaa  l.tll8 

probability  . Ffla  3.S0UB-02 


var 

design  point 

sensitivity  vectors 

x* 

U* 

alpha 

gasma 

delta 

eta 

sm 

4.S11E*05 

-3.076E-01 

-.1698 

-.1698 

.1722 

-.0590 

2.488E«01 

-5.378E-01 

-.2969 

-.2969 

.3098 

-.1800 

sw 

3 . 022E404 

8.876E-07 

.0000 

.0000 

.0000 

.0000 

SM 

4.9S9E406 

4.358B-01 

.2406 

.2406 

XU 

7.773E-01 

-1.484E«00 

-.8193 

-.8193 

.8193 

-1.2163 

xsw 

1 . 007E400 

1.332E-01 

.0735 

.0735 

-.0735 

-.0098 

xw 

9.920E-01 

6.818E-01 

.3763 

.3763 

-.3763 

-.2566 

xs 

1.155E-t>00 

1.496E-01 

.0826 

.0826 

-.0826 

-.0124 

»»  SENSITIVITY  ANALYSIS  AT  COMPONENT  LEVEL  «« 
type  of  parameters  for  sensitivity  analysis 


. . . iSVa  0 

isval  . distribution  parameters 

lsva2  . limit-state  fen  parameters 


isvaO  ..distribution  and  limit-state  fen  parameters 
sensitivity  with  respect  to  distribution  parameters 
limit -state  function  1 


d  (beta) /d  (parameter)  : 


var 

mean 

std  dev 

par  1 

par  2 

par  3 

par  4 

sm 

9.420E-06 

-3.22SE-06 

4.246E4^00 

-1.306E4^00 

tp 

1.709E-01 

-9.927E-02 

4.246E-f00 

-2.284E4^00 

sw 

-4.899E-07 

-4.349E-13 

-4.899E-07 

-4.349E-13 

mw 

-5.588E-07 

-1.201E-07 

0.OOOE4-OO 

0.000E4-00 

XU 

5.462E-»00 

-8.109E«00 

5.462E*00 

-8.109E^00 

xsw 

-1.471E+00 

-1.959E-01 

-1.471E+00 

-1.959E-01 

xw 

-2.788E>00 

-1.901E4-00 

-2.788E4>00 

-1.901E<»^00 

xs 

-2.394E«00 

-3.582E-01 

-2.394E>00 

-3.582E-01 

d(Pfl) /d(paraBMter)  : 

var 

mean 

std  dev 

par  1 

par  2 

par  3 

par  4 

sm 

-7.281E-07 

2.493E-07 

-3.282E-01 

1.009E-01 

fP 

-1.321E-02 

7.673E-03 

-3.282E-01 

1.765E-01 

sw 

3.787E-08 

3.361E-14 

3.787B-08 

3.361E-14 

mw 

4.319E-08 

9.283E-09 

O.OOOE-t-OO 

O.OOOE4.00 

XU 

-4.222E-01 

6.267E-01 

-4.222E-01 

6.267E-01 

xsw 

1.137E-01 

1.515E-02 

1.137E-01 

1.515E-02 

xw 

2.155E-01 

1.469E-01 

2.155E-01 

1.469E-01 

xs 

1.850E-01 

2.769E-02 

1.850E-01 

2.769E-02 

sensitivity  with  respect  to  limit-state  function  parameters 

limit -state  function  1 


par  d (beta) /d (parameter)  d(Pfl) /d(parameter) 

1  O.OOOE-fOO  0.000E4-00 


cyp«  oc  incagracxon  scmm  . . .  • 

itgal  . iaprovwd  Braltung  fonula 

ltga2  . laprovvd  BrAitung  fotaul* 

. k  Tv»dt'a  axacc  Intagral 

MX.  nuabcr  of  itoratlMS  for  ooch  fitting  point  ..inp«  4 


limit-stato  function  1 


coordinates  and  avs.  Min  curvaturas  of  fitting  points  in  rotated  space 


axis  u‘i 

u'n  G(u) 

u'i 

u'n  6(u) 

a'i 

1  1.810 

1.814  -4.0401-03 

-1.810 

1.814  -2.4161-03 

6.49931-04 

2  1.811 

1.812  -1.0171-04 

-1.811 

1.812  -8.2431-05 

1.09471-04 

3  1.812 

1.812  -2.9531-07 

-1.812 

1.812  -2.2011-07 

-3.79141-12 

4  1.812 

1.750  1.4911-04 

-1.812 

1.758  2.5541-04 

-1.75511-02 

S  1.812 

1.731  1.5141-04 

-1.812 

1.737  6.2081-04 

-2.36951-02 

6  1.811 

1.812  -1.4081-04 

-1.811 

1.812  -1.2961-04 

1.32901-04 

7  1.792 

1.831  -2.9681-01 

-1.790 

1.833  -1.7601-01 

6.29681-03 

improved  Breitung 

Tvedt's  11 

generalized 

reliability  index  betag  « 

1.7753 

1.7760 

probability 

Pf2  « 

3.7921-02 

3.7861-02 

Stop  -  Program  terminated. 
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